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ABSTRACT 

Although  simple  radial  equilibrium  turbine  design 
has  long  been  a  useful  design  method  for  axial  flow  gas 
turbine  blading,  more  demanding  requirements  for  high 
performance  machinery  make  such  an  approximation  less 
justified.   More  exact  analytical  methods  do  exist,  but 
require  many  times  more  effort,  even  if  a  computer  is 
utilized. 

A  relatively  simple  model  for  improving  a  simple 
radial  equilibrium  preliminary  design  is  proposed  in 
this  thesis.   The  axial  velocity  •distribution,  determined 
from  the  input  SRE  design,  is  modified  to  take  into 
account  the  results  of  flow  studies  in  flared  annuli 
and  of  actuator  disk  theory.   Using  the  corrected  axial 
velocity  distribution,  new  velocity  triangles  can  be 
calculated. 

A  FORTRAN  IV  computer  program  to  implement  the 
proposed  model  was  vs^ritten  and  is  presented.   Its 
potential  is  demonstrated  by  application  to  three 
example  simple  radial  equilibrium  designs.   The  results 
show  practical  differences  exist  between  a  simple  radial 
equilibrium  design  and  the  design  corrected  by  the 
proposed  method. 
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NOTATION 


A        flow  area,  ft 

a        coefficient  in  equations  2.24,  2.25,  2.27, 
and  2.2  8 

A        annulus  area,  ft^ 
ann 

b        coefficient  in  equations  2.24,  2.25,  2.27, 
and  2.2  8 

c        specific  heat  at  constant  pressure,  BTU/lbm/^R 

c        specific  heat  at  constant  volume,  BTU/lbm/°R 


g        32.2  ft  Ibm/lbf/sec 
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h  specific  enthalpy,  BUT/lbm;  actuator  disk 
height,  ft 

J  778  ft  Ibf/BTU 

M  Mach  number 

m  mass  flow  rate,  Ibm/sec 

N  turbine  rpm,  revolutions  per  minute 

P  pressure,  Ibf/ft^ 

Q  rate  of  heat  transfer  to  fluid,  BTU/sec 

r  radial  distance  from  machine  axis,  ft 

R  specific  gas  constant,  ft  lbf/lbm/°R 

s  specific  entropy,  BTU/lbm/°R 

T  temperature,  °R 

u  blade  velocity,  ft/sec 

V  fluid  velocity,  ft/sec 

• 

W  rate  of  output  shaft  work,  BTU/sec 

W  relative  velocity,  ft/sec 

X  axial  distance  from  an  actuator  disk,  ft 


a  absolute  angle  measured  from  axis,  degrees 

Aa  stator  turning  angle,  degrees 

3  relative  angle  measured  from  axis,  degrees 

AB  rotor  turning  angle,  degrees 

Y  specific  heat  ratio 

p  density,  Ibm/ft^ 

(J)  flare  angle,  degrees 

oj  angular  velocity,  radians/sec 

Subscripts 

h  hub  condition 

i  reference  condition 

o  total  (stagnation)  condition 

r  radial  component 

rel  relative 

t  tip  condition 

X  axial  component 

0  tangential  component 

1  rotor  inlet;  plane  one  condition 

2  rotor  outlet;  plane  two  condition 

3  plane  three  condition 

+0°  condition  far  downstream 

-°°  condition  far  upstream 

NOTE:   Notation  not  appearing  here  may  be  computer 
program  variables  listed  in  Appendix  A. 
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1.   INTRODUCTION 

1.1  General 

Energy  is  a  topic  of  increased  interest  in  current 
times.   Searches  for  new  sources  of  energy  and  new  means 
of  converting  energy  into  useful  forms  have  been  under- 
taken, along  with  efforts  to  conserve  the  resources 
available.   While  exploiting  solar  energy,  tapping 
geothermal  power,  and  drilling  for  oil  off-shore  receive 
much  public  interest,  progress  in  improving  existing 
conventional  forms  of  energy  transfer  and  conversion  must 
also  be  made.   The  great  importance  of  energy  utilization 
has  made  it  worthwhile  for  designers  of  energy  related 
systems  to  pay  even  m.ore  attention  to  achieving  the  best 
possible  performance.   Improving  the  efficiencies  of  the 
various  hardware  components  which  carry  out  the  heat 
transfers  and  energy  conversion  functions  of  the  thermal 
cycle  is  one  desirable  goal.   To  accomplish  this,  it  is 
useful  to  make  design  methods  more  accurate  and  more 
easily  applied  to  real  problems. 

1.2  Axial  Flow  Turbines 

This  thesis  deals  with  one  commonly  used  component 
for  converting  energy  into  useful  work--the  axial  flow 
turbine.   Turbines  convert  the  kinetic  energy  v;hich  has 
been  provided  the  v;orking  fluid  into  mechanical  work 
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performed  upon  a  rotating  shaft.   Axial  flow  turbines  are 
so  named  because  the  flow  of  the  \\7orking  fluid  through 
the  machine  is  predominantly  in  a  direction  parallel  to 
the  machine's  axis  of  rotation.   Such  turbines  can 
operate  with  a  variety  of  working  fluids,  and  in  either 
open  or  closed  thermal  cycles.   The  most  common  con- 
figurations are  as  a  steam  turbine  in  the  Rankine  Cycle 
(a  closed  cycle)  or  as  a  gas  turbine  in  the  Brayton 
Cycle  (an  open  or  closed  cycle) .   Historically,  the 
steam  turbine  found  earlier  application  than  the  gas 
turbine/  which  recently  has  gained  in  importance.   The 
development  of  the  gas  turbine  depended  upon  advances  in 
high  temperature  material  properties  and  was  also  closely 
tied  to  the  successful  operation  of  axial  compressors. 

1.3   Aerodynamic  Design  of  Axial  Flow  Turbines 
1.3.1   General 

One  aspect  of  the  design,  which  precedes  similarly 
for  axial  turbines  of  all  kinds,  is  the  aerodynamic 
design  of  the  turbine  blading.   The  efficiency  of  the 
turbine  depends  upon  the  smooth  transfer  of  kinetic 
energy  from  the  fluid  flow,  with  as  little  energy  wasted 
as  possible.   The  shapes  and  angles  of  the  turbine 
blading  must  be  appropriate  to  the  pattern  of  the  flow. 
The  aerodynamic  design  of  the  turbine  determines  what 
these  flows  and  shapes  should  be. 
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In  general,  the  aerodynamic  designer  of  an  axial 
turbine  follows  the  steps  listed  below: 

(1)  for  an  assumed  operating  point  and 
efficiency,  carry  out  the  thermodynamic  analysis  of  the 
thermal  cycle  to  determine  the  mass  flov;,  and  inlet  and 
outlet  states  for  the  turbine; 

(2)  postulate  a  working  fluid  flow  which 
satisfies  the  laws  of  motion  and  produces  the  required 
work  ; 

(3)  design  the  shape  of  the  turbine  blading 
appropriately; 

(4)  calculate  the  efficiency  and  evaluate  the 
design; 

(5)  repeat  all  of  the  above  steps  as  necessary. 
1.3.2   Flow  Solutions 

As  is  true  of  most  design  efforts,  there  are  many 
levels  of  sophistication  and  detail  which  the  designer 
may  seek,  depending  upon  the  stage  of  the  design  and/or 
the  cost  and  criticality  of  the  product.   In  the  aero- 
dynamic design  of  a  turbine,  the  designer  must  decide 
which  level  is  best  in  each  case.   For,  although  the 
equations  of  motion  and  continuity  for  a  fluid  are 
readily  v/ritten,  they  cannot  be  solved  in  practice  unless 
simplifying  assumptions  are  made.   Fortunately,  such 
simplified  forms  have  proved  adequate  in  many  applications, 
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1.3.3  Two -Dimensional  Flov^? 

If  the  hub-tip  ratio  of  a  turbine  is  nearly  1.0, 
no  large  radial  component  of  fluid  velocity  is  expected, 
and  flow  conditions  will  be  approximately  the  same  from 
root  to  tip.   The  flow  can  then  be  analyzed  in  two 
dimensions.   This  results  in  blades  of  constant  cross- 
section,  v;hich  have  been  extensively  studied  experimentally 
in  the  form  of  cascades. 

1.3.4  Three-Dimensional  Flow 

In  some  cases  the  designer  faces  requirements  on 
blade  tip  speed,  pressure  ratio,  and  mass  flow  which 
could  result  in  low  hub-tip  ratios  (i.e.,  relatively 
long  blades)  and  considerable  turbine  flare.   These 
conditions  reduce  the  accuracy  of  simple  two-dimensional 
design  methods  and  force  the  consideration  of  three- 
dimensional  flov7.   The  general  solution  is  impractical, 
but  simplified  methods  have  been  developed  and  used  with 
success.   Several  of  these  approaches  are  mentioned  below. 

1.4   Approaches  to  Three-Dimensional  Design 

The  first  assumption  made  in  virtually  all  cases  is 
to  assume  axial  symmetry.   The  number  of  blades  is 
sufficiently  large  that  circumferential  conditions  at  a 
station  can  be  represented  by  average  values. 
1.4.1   Radial  Equilibrium 

Assuming  axisymmetric  flow,  the  differential 
equation  of  motion  in  the  radial  direction  for  the  fluid 
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can  be  put  in  a  form  (see  paragraph  2.4.1)  which  may  be 
numerically  integrated  by  a  computer.   The  procedure  in- 
volves complicated  iterations,  however,  and  elaborate 
measures  must  be  taken  to  ensure  the  proper  convergence 
of  the  streamline  locations.   Carmichael  has  discussed 
the  procedure  and  presented  block  diagrams  of  possible 
computer  programs.    This  method  allovv'S  for  fluid 
velocities  in  the  radial  direction  which  result  from 
imbalances  in  the  static  pressure  and  centrifugal 
acceleration  forces  on  the  fluid,  but  its  complexity 
is  formidable. 

1.4.2  Simple  Radial  Equilibrium 

If  the  radial  components  of  the  fluid  velocity  are 
neglected,  the  so-called  streamline  curvature  term  in  the 
differential  equation  drops  out.   The  differential 
equation  of  motion  in  the  radial  direction  can  then  be 
written  involving  only  the  axial  and  tangential  velocities 
This  can  be  easily  integrated  for  axial  velocity  when 
given  a  reasonable  distribution  for  the  tangential 
velocity  (see. paragraph  2.4.3).   Simple  pairs  of  functions 
.for  axial  and  tangential  velocity,  so  determined,  have 
been  widely  used  in  three-dimensional  design. 

1.4.3  Actuator  Disk  Theory 

Simple  radial  equilibrium  analysis  depended  upon  the 
radial  fluid  velocity  being  zero  before  and  after  each 
blade  row.   In  most  cases,  changes  in  the  radial  distri- 
bution of  density  result  in  radial  shifts  in  the 
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streamlines.   Therefore,  a  non-zero  radial  velocity 
component  must  exist. 

An  alternative  to  the  simple  radial  equilibrium 
analysis  is  the  actuator  disk  formulation.   Each  blade 
row  is  replaced  by  a  narrow  "actuator  disk"  which  provides 
the  appropriate  sudden  change  in  fluid  tangential  velocity 
This  device  permits  an  approximate  solution  to  the 
equations  of  motion,  allowing  finite  radial  velocities, 
but  still  is  restricted  to  annular  flow  between  walls 
parallel  to  the  machine  axis.   (See  paragraph  2.5) 

1.5  Fluid  Flow  in  Flared  Annuli 

A  frequent  design  practice  is  to  keep  the  axial 
velocity  approximately  constant  throughout  the  machine, 
compensating  for  changes  in  fluid  density  by  increasing 
the  annulus  area.   This  area  change  results  in  "flare"; 
that  is,  the  hub  and  shroud  surfaces  are  not  parallel  to 
the  machine  axis,  but  are  sloped. 

This  flare  could  be  expected  to  cause  a  general 
change  in  the  axial  velocity  distribution  within  the 

annulus,  with  significant  radial  velocity  components 

2 

appearing  near  the  walls.   Carmichael  and  Par   have 

suggested  a  relation  which  can  be  used  to  calculate  the 
change  in  axial  velocity  in  flared  annuli. 

1.6  Summary 

The  preceding  discussion  was  intended  to  demonstrate 
a  heed  for  a  relatively  simple  method  to  predict  the 
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fluid  flow  distribution  in  the  presence  of  flare,  improving 
upon  the  accuracy  of  the  simple  radial  equilibrium  model. 
Such  a  method  would  hopefully  be  useable  in  the  preliminary 
design  phase  with  little  increased  effort,  as  compared  to 
the  more  complicated  streamline  curvature  approaches. 

A  method  is  proposed  in  this  thesis  for  "correcting" 
a  preliminary  simple  radial  equilibrium  turbine  design 
for  flare  and  strefimline  curvature  effects.   This  is 
accomplished  by  modifying  the  axial  velocity  radial 
variation,  as  determined  from  the  simple  radial  equilibrium 
equations.   A  correction  for  turbine  flare  is  applied 
using  the  findings  of  Pai;  next,  the  axial  velocity  is 
further  modified  at  the  blade  leading  and  trailing  edges 
by  the  approximate  solutions  to  the  actuator  disk  analysis. 
As  the  computations  are  still  iterative  in  nature,  they 
are  very  suitable  for  computerization. 

A  computer  program  to  modify  a  sim.ple  radial  equili- 
brium turbine  design,  according  to  the  model  described 
above,  is  presented  in  this  thesis.   While  not  theoreti- 
cally deep-seated,  it  is  hoped  the  methodology  may  be 
helpful  in  some  turbine  design  applications. 


17 


2 .   BACKGROUND 

The  purpose  of  this  chapter  is  to  provide  the  back- 
ground for  the  methodology  employed  in  the  computer  model 
to  be  discussed  in  Chapter  3.   The  emphasis  is  placed  on 
those  assumptions  and  relations  which  directly  apply  to 
the  program  model,  although  some  additional  development 
is  included  for  perspective  and  to  point  to  possible 
areas  of  application  and  further  work. 

2.1   Fluid  Mechanics 

Axial  turbine  flow  is  frequently  modeled  as  incom- 
pressible, at  least  in  the  radial  direction,  so  that  the 
equations  of  motion  for  an  incompressible  fluid  apply. 
It  is  also  known  that  viscous  effects  can  usually  be 
neglected  outside  the  boundary  layer.    In  cylindrical 

coordinates,  the  equations  of  motion  for  an  ideal,  incom- 

3 

pressible  fluid  can  be  written  as 

V  _£  +  -A^  +v-^--^=--^|^  (2.1) 

r  dr  rd9      x  9x     r       p     6r 

r  8r    r  89      x   9x     r      pr  dd 

3V     V„  8V        9V       g   .„ 

V_2i  +  -l-^4-V   -3ii=-^|i,  (2.3) 

r  3r    r   96     x   9x      p  9x 
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where  V  ,  V^ ,  and  V^  refer  to  components  of  fluid  velocity 
in  the  radial,  tangential,  and  axial  directions,  respect- 
ively.  P  is  the  static  pressure.   The  continuity 
equation  is 


8V     V      ,   9V„    8V 

r  ,   r   ,  1   0  ,    X    _  /o  /1^ 


2.2   Thermodynamics 
2.2.1   Work 

From  the  first  law  of  thermodynamics,  the  steady 
flow  energy  equation,  describing  the  change  in  state  of 
a  fluid  system  in  flov;ing  through  a  control  volume,  can 
be  written  as 


Q  -  W  =  m[(h  +  ^ -)   ^  -  (h  +  "^ -)  .  ]  ,         (2.5) 

2g  J  out        2g  J  m 
^o  ^o 


where  potential  terms  have  been  ignored.   Q  is  the  rate 
of  heat  transfer  to  the  fluid,  W  is  the  rate  of  output 
shaft  work,  m  is  the  mass  flow,  V  is  the  fluid  velocity, 
and  h  is  the  enthalpy.   A  useful  definition  is  to  make 
stagnation  enthalpy,  h  , 


h   =  h  +  ^^  (2.6) 

o        2g  J 


Using  equation  (2.6),  and  assuming  that  the  heat  transfer, 
Q,  from  an  axial  turbine  is  small  (compared  to  W) , 
equation  (2.5)  becomes 
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-«=  ^\Z    -    ^ol'-™'  (2.7) 


where  the  2  and  1  refer  to  the  outlet  and  inlet  states 
respectively. 

Following  reference  1,  the  rate  of  output  shaft  work, 
W,  is  also  equal  to  the  product  of  the  angular  velocity, 
CO,  of  the  blades  and  the  torque  caused  by  the  moments  of 
the  external  reaction  forces.   The  torque  is  given  by  the 
change  in  angular  momentum  of  the  fluid,  or 


Torque  =  mir^V^^    ~    ^i^ei^ *  (2.8) 


Since 


W  =  Torque  •  co 


and  the  blade  velocity,  u,  is  just  ojr ,  then 


•g  J 
m       ^o 


(2.9) 


In  axial  machines  it  is  a  good  assumption  to  make  U2=u,=u 
for  preliminary  design. 
2.2.2   Gas  Lav/s 

The  equation  of  state  for  a  perfect  gas  can  be 
written 
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P  =  PRT,  (2.10) 

where  R  is  the  specific  gas  constant.   The  following 
relations  between  the  specific  heats  are  the  standard 
ones. 


C 

'  -  C 

V 


C   -  C   =  R/J 
p     V 


Cp  =  (^)  §.  ■  (2.11) 


2.2.3   The  Second  Law 

From  the  second  law  of  thermodynamics,  using  Gibbs 
Equation  and  the  perfect  gas  law,  the  specific  entropy, 
s,  can  be  written  as  ' 


Tds  =  dh  -  -^  dP.  (2.12) 

pu 


2.3   Compressible  Flow  Relations 

In  all  but  very  low  speed  turbomachines ,  in  which 
velocity,  pressure,  and  temperature  changes  are  very 
small,  the  effects  of  compressibility  of  the  working 
fluid  become  important.   Consideration  of  three-dimen- 
sional com.pressible  flow  is  extremely  complex.   Fortunately, 
the  changes  in  fluid  density  are  often  not  great  in  the 
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radial  direction  at  a  station,  and  need  to  be  considered 
only  in  the  axial  flov;  direction. 

The  one-dimensional  compressible  flow  relations  for 

an  ideal  gas  appear  in  many  representations  and  are 

3 

conveniently  used.   Horlock   has  a  discussion  of  their 

derivation  and  summarizes  the  results.   The  following 
forms  will  be  useful  in  discussion  later.   The  meanings 
of  the  symbols  can  be  found  in  the  notation  list.   From 
equation  (2.6),  since,  for  a  perfect  gas,  h  =  c  'T, 


T   =  T  +  J^'    T^  .  (2.13) 

o       2g  JC 
o   p 


Using  (2.11)  and  dividing  by  T  ,  equation  (2.13)  can  be 
rewritten  as 


1=1-  (^)— ^.  (2-14) 

■^O  '   2g  R  T 

^o    o 


The  Mach  number,  M,  is  defined  by 


M  =  -  =  ^-^ .  (2.15) 

^    / ~ 

g  YRT 


From  equations  (2.14)  and  (2.15), 


,       T  g  yR      t 
(1)2  =   _o^o_[_  _  xzl,  (2.16) 
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The  pressure-temperature  relation  for  an  isentropic 
process  can  be  shown  to  apply  in  the  compressible  case, 
and  gives 


^=  (^)^-^  (2.17) 

o     o 


From  continuity,  the  flow  area.  A,  is 


A  =  ^,  (2.18) 


which  can  be  determined  from  equations  (2.14)  and  (2.17) 
above . 

2.4   Simple  Radial  Equilibrium 
2.4.1   General 

As  indicated  in  the  introduction,  when  the  hub-tip 
ratio  of  a  turbine  stage  is  further  and  further  from  the 
value  of  1.0  (i.e.,  relatively  longer  and  longer  blades), 
the  two-dimensional  flov/  design  methods  become  less  and 
less  adequate.   In  some  marginal  cases,  three-dimensional 
details  of  the  design  are  ignored  in  favor  of  other 
considerations,  for  example,  maintaining  constant  blade 
sections.   But  modern  performance  requirements  make  such 
practices  less  attractive.   The  general  three-dimensional 
flow  equations  are  difficult  to  solve,  however,  and 
simplifying  assumptions  are  frequently  made.   The  most 
general  such  assumption  is  to  neglect  any  variation  of 
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properties  in  the  circuraferential  direction.   In  the 
axial  direction,  where  fluid  property  changes  may  be 
large,  the  one-dimensional  compressible  flow  relations 
are  assumed  to  describe  the  flow.   Finally,  in  the  radial 
direction  the  equation  of  motion  in  a  simplified  form-- 
termed  the  simple  radial  equilibrium  (SRE)  equation--is 
often  applied.   A  discussion  of  the  development  and  use 
of  the  simple  radial  equilibrium  equation  follows. 
2.4.2   The  Simple  Radial  Equilibrium  Equation 

Assuming  axial  symmetry,  the  equation  of  motion  in 
the  radial  direction,  equation  (2.1),  can  be  written  as 


8V        dV  V  2      g 

V   -^+V   _^--i-=--^|£.  (2.19) 

r   8r     X   8x     r       p   3r 


From  the  identity  V^=V^+V.^+V^  and  by  differen- 

•^       r     e     X       -^ 

tiating  with  respect  to  r,  equation  (2.6)  becomes 


.,    3h      ,       dV  8V       9V 

^=^--^[V   _^  +  v„.^+V   3-^],      (2.20) 
dr        dr  q  J   r  9r     G    3r    x  dr 

o 


and  equation  (2.13)  can  be  rewritten  in  the  form 


1   |P  =  t|^-|11.  (2.21) 

pj  9r      3r    9r 


Substitution  of  equations  (2.20)  and  (2.21)  into  equation 
(2.19)  gives 

.^  dh  8V       9V       9V_    V„^ 

t|^  — ^  =  V  —I.  -   V   -^   -    V^-^l  -  -iL..  (2.22) 

9r    9r     x  9x     x  9r     6  9r     r 
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At  this  point  9s/9r  and  3h  /3r  are  often  assumed  to  be 

zero.   According  to  Horlock   these  two  assumptions, 

which  prove  to  be  nearly  true  in  practice,  hold  for  the 

flow  between  blade  rows  of  a  reversible  turbine  in  which 

delivered  work  and  total  pressure  drop  across  the  blade 

are  the  same  at  all  radii.   The  term,  V  • 9V  /8x,  is 

X    r 

called  the  streamline  curvature  term  and  is  the  cause  of 
the  difficulty,  mentioned  in  the  introduction,  in 
obtaining  a  solution  to  equation  (2.22).    In  addition 
to  assuming  3s/3r  and  3h  /3r  are  equal  to  zero,  if  it  is 
assumed  that  V   equals  zero  at  the  leading  and  trailing 
edges  of  a  blade  row,  the  streamline  curvature  term  drops 
out  and  the  partial  derivatives  can  be  replaced  by 
ordinary  derivatives,  giving 


dV^         dV,    V  2 

V  —21  =  -  V  •  -^  -  -^.  (2.23) 

X   dr       9    dr     r 


Equation  (2.23)  is  called  the  simple  radial  equilibrium 
equation. 

A  more  direct  derivation  of  the  simple  radial 
equilibrium  equation,  which  may  be  more  intuitive,  can  be 
developed  from  consideration  of  an  infinitesimal  volume 
of  fluid  in  equilibrium  with  the  external  forces  acting 
on  it.    The  details  are  shown  in  Figure  2-1,  on  the 
next  page.   The  substitution  for  the  g  /p(dP/dr)  term 
precedes  as  above. 
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dr 


A'W 


\  / 

V 


Sum  of  forces  in  the  radial  direction  =  0 

dP  ^ 

r.d9(P  +  ~  dr)  -  P-r-dS-  — (pr-d9-dr)  =  0 


riP  ft 

~  r-de-dr  =p-|  r-dG-dr 


1  ^  =  _^ 
p  dr  ~     r 


Figure  2-1 
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In  summary,  the  simple  radial  equilibrium  equation 
is  valid  for  the  following  conditions: 

(1)  ideal,  incompressible  fluid  flow 

(2)  axial  symmetry 

(3)  perfect  gas  law  holds 

(4)  V  =0  before  and  after  blade  rows 

r 

(5)  dh  /dr  =0 

(6)  ds/dr  =  0 

2.4.3   Solutions  to  the  Simple  Radial  Equilibrium  Equation 

From  equation  (2.9)  it  can  be  seen  that  the  work 
performed  by  an  axial  turbine  is  related  to  the  tangen- 
tial velocities  before  and  after  the  rotor,  V„,  and  V^^, 
respectively.   A  design  problem  for  an  axial  turbine  is 
to  propose  a  reasonable  tangential  velocity  distribution 
which  will  give  the  desired  work,  and  then  by  means  of 
the  simple  radial  equilibrium  equation  find  the  axial 
velocity,  V  ,  so  that  the  blade  shapes  can  be  determined. 
For  simple  functions  of  V^^  =  V^  (r)  ,  equation  (2.23)  is 
easily  integrated  directly  for  V  =  V^(r). 

It  would  appear  that  independent  functions  for  V„, 
and  V„2  could  be  chosen  and  corresponding  functions 
V  , (r)  and  V  ^ (r)  calculated  from  the  simple  radial 
equilibrium  equation.   However,  the  V^ ' s  are  related  by 
the  v/ork  relation,  equation  (2.9),  and  by  the  assumption, 
used  in  deriving  the  simple  radial  equilibrium  equation, 
that  the  stagnation  enthalpy,  h  ,  after  each  blade  row, 
does  not  vary  with  radius. 
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A  familiar,  useful  set  of  tangential  velocities 
given  by  Carmichael  and  Lewis,  and  described  by  Horlock  , 
are  the  following: 


V0j_  =  ar"  +  I  (2.24) 


Vg2  =  ar^  -  I  (2.25) 


where  a,  n,  and  b  (b  >  0)  may  be  chosen  by  the  designers. 
Then,  from  equation  (2.9) 


W  =  E_  (V    -  V   )  =  ^^^^   (V    -  V   ) 

g  J  ^^61    ^02^    60g  J  ^^61    ^62^ ' 
m    ^o  ^o 


and   using   equations    (2.24)    and    (2.25) 


W         27TrN       ,2b,  47TNb  ,^    ^^. 

m  -^o  ^o 


Equation  (2.26)  shows  the  work  done  (hence  h  )  is 
independent  of  radius  and  has  the  proper  sign  for  a 
turbine . 

A  modification  of  the  set  of  tangential  velocities 
(2.24)  and  (2.25)  is  considered  in  this  thesis.   Suppose 
a  set  of  tangential  velocities  of  the  following  form  is 
given: 
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^1 

^61  "  ^2^   "*"  ^1^  "^  ^o  "^  r  '  (2.27) 

b 


Here  also  the  work  is  independent  of  radius,  for 


m     ^o 


Here  b,  must  be  greater  than  b„  if  work  is  to  be  done  by 
the  turbine.   The  constant  b,  is  virtually  always  positive, 
and  b„  is  usually  negative  or  zero.   The  advantage  of 
allov/ing  b,  and  b^  to  be  different  is  to  be  able  to 
provide  for  special  cases  such  as  -free  vortex,  axial 
outflow  (a  ,  a,  ,  a  ,  b_  =  0;  b,  7^  0)  . 

Substituting  an  expression  of  the  form  of  (2.27) 
into  the  simple  radial  equilibrium  equation  and  integrating 
from  reference  point  V  .  =  V  (r.)  to  point  V  =  V  (r) ,  the 

^         XI      X    1        '^  XX 

following  expression  for  axial  velocity  is  obtained  (the 
details  are  in  Appendix  D) : 


(4a^a2  +  2a^2  j  (^2_^  _^2  ^  .g  (^^^^^^^b^)  (r-r^) 


(2a  2+4a,b,  )  In  (— )  +  2a  b,  (-  -  — )  .         (2.30) 
^   o    11    r .      o  1  r   r . 

1  1 
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Equations  (2.27)  and  (2.30)  for  tangential  velocity  and 
axial  velocity  form  a  set  which  satisfy  the  simple  radial 
equilibrium  equation  (2.23).   It  is  important  to  note 
that  many  factors  and  various  criteria  for  good  practice 
influence  the  detailed  choice  of  the  constants  in  the 
expressions  for  V^. 

Tv;o  well  knov/n  special  cases  of  equations  (2.27) 
and  (2.28)  are  the  constant  reaction  V^   distribution  and 
the  exponential  Vq  distribution.   These  tangential 
velocity  distributions  are  presented  in  Appendix  D 
together  with  the  corresponding  function  for  axial 
velocity. 

2.5   Actuator  Disk  Theory 
2.5.1   General 

It  is  recalled  that  the  basis  for  the  simple  radial 
equilibrium  design  method  is  that  the  radial  component  of 
the  fluid  velocity  be  zero  at  the  leading  and  trailing 
edges  of  a  blade  row.   This  assumption  results  in  stream- 
lines of  the  form  illustrated  in  Figure  2-2  (a) .   A  more 
general  method  of  analysis,  which  takes  some  account  of 
radial  shifts  in  the  streamlines  between  blades,  is  the 
actuator  disk  model.   As  described  by  Horlock  ,    a  rotor 
or  stator  rov/  is  replaced,  conceptually,  at  approximately 
midchord,  by  a  thin  disk  across  which  the  tangential 
velocity  change  takes  place.   Such  a  device  allows  the 
equations  of  motion,  including  radial  velocity  terms,  to 
be  solved  analytically,  although  the  formulation  is  still 

complex. 
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2.5.2   Approximate  Solution 

With  some  simplifications,  an  approximation  to  the 
complete  actuator  disk  results  can  be  obtained.   Radial 
velocities  are  assumed  to  be  zero  at  the  walls  and  small 
elsewhere,  and  simple  radial  equilibrium  is  assumed  to 
be  valid  a  sufficient  distance  from  the  disk.   According 
to  this  model,  the  streamlines  will  be  similar  to  those 
shown  in  Figure  2-2  (b) .   The  approximate  solutions  for 
the  axial  velocity  near  a  single  disk,  in  terms  of  the 
axial  velocity  on  the  same  streamline  far  up  and  down- 
stream, are: 


(V  )  ,   -  (V  ) 
V^(x)  =(V^)_^  +  -^^^ exp  C-^),        (2.31) 


upstream  of    the    disk    (x   <    0),    and 


(V    )         -     (V    ) 
^x'-'    =    <\'+~  -   ^^^ ^^  exp(-I^),  (2.32) 


downstream  of  the  disk  (x  >  0) .   The  notation  (V  ) 

and  (V  )    refer  to  the  axial  velocities  given  by  simple 
X  —  °°    ' 

radial  equilibrium  a  long  way  downstream  and  upstream, 
respectively;  h  is  the  height  of  the  actuator  disk. 

This  analysis,  like  the  simple  radial  equilibrium 
analysis  described  in  the  previous  section,  treats  ideal, 
incompressible  flow  in  cylindrical  annuli. 
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For  real  turbines,  where  many  rows  of  blades  (disks) 
may  influence  the  axial  velocity  distribution,  the  effects 
of  each  disk  may  be  considered  separately,  then  super- 
imposed to  find  the  real  axial  velocity.    This  procedure 
is  used  in  Appendix  E  to  find  the  axial  velocity  distri- 
bution at  a  point  influenced  by  two  nearest  adjacent 
disks. 

2.6   Flow  in  Flared  Annuli 
2.6.1   General 

In  addition  to  improving  the  methods  for  cylindrical 
flow,  efforts  to  understand  the  fluid  flow  in  ducts  of 
changing  annular  cross-section  have  been  made.   Such 
sections  are  commonly  used  in  practice. 

Korlock,  in  reference  6,  describes  some  theoretical 
models,  proposed  by  various  workers,  for  the  determination 
of  axial  velocity.   Analyses  by  Lewis  and  Horlock  for 
incompressible  flow  in  conical  ducts,  by  Wu  for  compressible 
flov;  in  flared  annuli,  and  a  simplification  of  Wu's  analysis 
by  Walker  are  discussed.   The  drawback  of  these  analytical 

methods  is  their  complexity. 

2 

Experimental  work  by  Carmichael  and  Pai   has 

suggested  a  simple  expression  relating  the  axial  velocity 
of  flow  in  a  flared  duct  to  the  angle  the  streamline 
makes  with  the  axis.   Their  work,  described  below,  forms 
the  basis  of  the  flare  correction  to  axial  velocity 
given  in  this  thesis. 
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2.6.2   Experiments  Representing  Fluid  Flow  in  Flared 

2 

Turbines , 

Experiments  were  conducted  using  a  tilted-floor 
electrolytic  tank  to  measure  flows  within  a  flared 
annulus.   An  electrolytic  tank  provides  an  analogue  to 
ideal  incompressible  fluid  flow.   Electrical  potential 
measured  in  the  tank  is  analogous  to  fluid  potential  and 
the  electrical  potential  gradient  is  analogous  to  the 
fluid  velocity.   A  sector  of  an  axisymmetric  model  can 
be  represented  by  tilting  the  floor  of  the  tank. 

Although  flow  in  a  real  turbine  is  compressible, 
the  radial  variation  in  density  at  any  axial  station  is 
relatively  small.   Consequently,  it  can  be  assumed  that 
in  usual  geometries  the  strearaline  positions  for  the 
real,  compressible  fluid  are  very  nearly  the  same  as  for 
the  ideal,  incompressible  case. 

One  aspect  of  the  fluid  flow  measured  was  the 
general  change  in  axial  velocity  distributions  in  the 
flared  annulus.   Three  models  were  studied,  with 
different  distributions  of  flare  between  the  hub  and 
shroud  walls;  each  model  had  a  hub-tip  ratio  of  0.9. 
The  distribution  was  expected  to  be  a  function  of  location 
within  the  annulus,  described  by  an  angle,  <{) ,  defined  as 
shown  in  Figure  2-3.   The  flow  conditions  were  measured 
at  an  axial  position,  free  of  discontinuities,  such  as 
section  A-A'  in  Figure  2-4.   The  axial  velocities  were 
found  to  have  the  distribution  shown  in  Figure  2-5, 
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approaching  the  flow  from  a  source.   The  equation, 
according  to  potential  flow  theory,  for  the  axial 
velocity  distribution  for  a  two-dimensional  source  is 


V  . 


Pai  suggests  that  comparison  of  this  equation,  also 
plotted  in  Figure  2-5,  with  the  experimental  results, 
indicates  that  for  angles  of  (})  <  ±  15°,  equation  (2.33) 
is  a  good  representation  of  the  axial  flow.   Even  at 
larger  angles,  the  agreement  may  be  good  enough  to  use 
equation  (2.33)  for  some  design  applications. 

2.7   Axial  Turbine  Flow  Geometry 
2.7.1   Velocity  Triangles 

It  should  be  restated  that  the  purpose  of  deter- 
mining the  details  of  the  fluid  flow  is  to  be  able  to 
design  the  blades  appropriately.   The  relevant  fluid 
velocities  and  flow  angles  are  conveniently  displayed 
in  the  form  of  vector  diagrams,  often  called  velocity 
triangles.   Velocity  triangles  for  a  single  stage  axial 
turbine  with  pure  axial  inflow  are  shown  in  Figure  2-6, 
where  any  radial  velocity  components  have  been  neglected 
The  calculation  of  such  velocity  triangles  for  each 
leading  and  trailing  edge  and  at  several  radii  must  be 
accomplished  before  the  detailed  design  of  the  turbine 
blading  can  continue. 


37 


Stator  inlet 


Rotor  outlet 


Rotor  inlet 


H— Ve2 


Remarks : 


1,  the  signs  of  velocit:/  corrponents  are  positive  in  the 
direction  of  U  or  V^ 

2,  angles  are  positive  in  a  coimtercloclc.d.se  direction 
from  V„ 


3.     as  dravm  above — 


V0-| ,  OC^ ,    /S.    are  positive 
Vq2,  (X2j     B>2  ^^®  negative 


Figure  2-6 
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2.7.2   Annulus  Area 

Suppose  a  fluid  is  flowing  uniformly  through  a  small 
area  at  a  given  rate.   If  the  area  encompassing  the  flow 
is  perpendicular  to  the  streamlines,  it  is  termed  the 
flow  area,  A,  and  is  given  by  equation  (2.18).   If  the 
area  is  in  a  surface  whose  perpendicular  makes  an  angle 
a  to  the  streamlines,  as  in  Figure  2-7,  then  the  area 
required  to  encompass  the  same  flow  is 


A    =  A/cos  a.  (2. 34) 

ann 


In  the  case  of  an  axial  flow  turbine,  A    represents 

ann 

the  annulus  flow  area  perpendicular  to  the  machine  axis, 
where  a  is  the  same  a  as  in  Figure  2-6.  From  Figure  2-7 
it  can  be  seen  that 


cos  a   =  V  /V; 


whence  equation  (2.18)  becomes 


m 


ann     pV 

X 


(2.35) 
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3.   MODEL  DESCRIPTION  AND  IMPLEMENTATION 

3.1  General 

In  the  introductory  chapter  the  need  for  a  simple 
methodology  to  generate  a  preliminary  axial  turbine  design, 
taking  account  of  turbine  flare  and  radial  streamline 
shifts,  v;as  discussed.   The  simple  radial  equilibrium 
design  method  was  presented  in  Chapter  Two,  together 
with  the  results  of  actuator  disk  theory  and  axial  flow 
distribution  studies  in  flared  annuli.   In  this  thesis 
a  procedure  is  proposed  which  accepts  a  simple  radial 
equilibrium  preliminary  design  as  a  basis,  and  then 
modifies  the  axial  velocity  distribution  in  accordance 
with  the  actuator  disk  and  flare  corrections  mentioned 
above.   These  modifications  to  the  simple  radial  equili- 
brium flov;  solution  result  in  slightly  different  blade 
shapes,  which  should  be  more  commensurate  with  the  real 
flow. 

This  chapter  will  present  the  model  for  modifying  a 
simple  radial  equilibrium  preliminary  design,  and  describe 
a  computer  program  which  implements  the  model. 

3.2  Model  Description 

To  begin,  it  is  useful  to  consider  a  sketch  of  the 
geometry  and  then  a  bare  outline  of  the  problem. 


41 


Figure  3-1  represents  a  single  stage  axial  flow  turbine. 
The  stator  inlet,  rotor  inlet,  and  rotor  outlet  planes 
are  labeled  1,  2,  and  3  respectively.   Other  features 
will  be  discussed  as  they  become  important. 

The  following,  are  steps  in  the  procedure  which  may 
be  identified  in  order  to  guide  the  discussion: 

(1)  "input"  a  simple  radial  equilibrium 
preliminary  design 

(2)  at  each  of  the  three  principal  planes  of  the 
turbine  stage,  calculate  a  set  of  streamlines  which  is 
consistent  v.'ith  certain  given  conditions  and  the 
modified  axial  velocity 

(3)  at  every  streamline  in  each  principal 
plane,  calculate  velocity  triangle  data. 

The  rest  of  this  section  will  discuss  the  above  steps  in 
general  terms.   The  details  of  the  procedure  are  best 
understood  in  the  context  of  the  computer  program  des- 
cription, beginning  in  paragraph  3.3. 
3.2.1   The  Input  Design 

The  purpose  of  the  procedure  discussed  in  this 
thesis  is  to  make  slight  "corrections"  in  an  existing 
simple  radial  equilibrium  design  in  order  to  refine  the 
detailed  blade  shapes.   As  a  result  of  carrying  out  a 
preliminary  design  for  a  single  stage,  many  parameters 
will  have  been  established.   The  working  fluid,  total 
temperatures  and  pressures,  hub  and  tip  radii,  machine 
rpm,  mass  flov/,  and  blade  widths  are  known;  also,  in  the 
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case  of  simple  radial  equilibrium,  an  appropriate  pair  of 
functions  for  axial  and  tangential  velocity  at  each  plane 
will  have  been  picked.  These  quantities  are  indicated  in 
Figure  3-1. 

For  present  purposes,  it  is  assumed  that  it  is 
desired  to  retain  the  input  annular  geometry  (hub  and 
tip  radii,  and  blade  widths)  in  the  modified  case, 
although  other  choices  are  possible.   For  example,  the 
average  axial  velocity  could  be  fixed,  and  the  tip  radius 
allowed  to  change  to  accorainodate  the  flow. 

For  simplicity,  the  input  preliminary  design  is 
restricted  to  having  hub  and  shroud  surfaces  of  uniform 
slope  over  the  v;hole  stage,  as  in  Figure  3-1.   Also,  the 
rotor  tip  clearance  and  stator-rotor  gap  are  assumed  to 
be  zero.   The  model  could  be  applied  to  more  general 
cases  v;ith  relatively  minor  modifications  to  the  computer 
program.   The  turbine  flare  may  be  distributed  between 
the  hub  surface  and  shroud  surface  in  any  proportion. 

It  should  be  noted  that  the  procedure  described 
here  will,  in  no  sense,  improve  a  bad  preliminary  design. 
3.2.2   Calculating  the  Streamlines 

To  fully  determine  the  blade  shape,  fluid  flow 
information,  including  axial  velocity,  must  be  known  at 
several  radii.   If  the  design  were  merely  simple  radial 
equilibrium,  the  velocity  triangles  could  be  independently 

found  at  any  convenient  radius,  since  V   and  V^  (and  U) 

-^  x      y 

are  knov/n  functions  of  r  alone.   However,  in  the  present 
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case,  the  actuator  disk  correction  to  the  axial  velocity 
at  a  location  on  one  plane  is  dependent  upon  the  axial 
velocity  along  the  same  streamline  at  stations  up  and 
downstream.   For  this  reason  a  set  of  streamlines  passing 
all  three  principal  planes  must  be  calculated. 

To  find  the  final  set  of  streamline  positions,  three 
different  iterative  loops  are  involved.   They  are: 

(1)  at  a  given  plane  and  streamline,  iteration 
is  required  to  obtain  the  next  (radially  outward) 
streamline  position; 

(2)  at  a  given  plane,  iteration  is  required  to 
obtain  an  average  axial  velocity  which  will  result  in  a 
match  between  the  fixed  mass  flow  rate  and  annulus  area; 

(3)  since  a  change  in  the  axial  velocity  at  a 
streamline  in  one  plane  affects  the  axial  velocity  in 
another  plane,  the  streamlines  are  recalculated  at  each 
plane  in  succession  until  the  changes  from  one  iteration 
to  the  next  become  insignificant. 

The  modified  axial  velocity,  used  in  the  streamline 
calculations,  is  found  utilizing  the  approximate  actuator 
disk  results  described  in  paragraph  2.5  and  Appendix  E, 
and  the  flared  annulus  axial  flow  relation  described  in 
paragraph  2.6.   The  only  changes  to  the  simple  radial 
equilibrium  preliminary  design  are  a  result  of  these  two 
"corrections"  to  the  axial  velocity,  which  impact  on  the 
velocity  triangles.   The  changes  represent  a  departure 
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from  simple  radial  equilibrium.   The  justification  for 
this  procedure  is  that  the  adjustments  are  small  and  are 
improvements  in  the  direction  of  the  real  solution. 
3.2.3   Velocity  Triangle  Data 

After  the  positions  of  streamlines  have  converged, 
all  relevant  parameters  in  the  velocity  triangles  can  be 
calculated  from  the  primary  variables  U,  corrected  V  , 
and  Vq ,  and  input  parameters.   Other  useful  design  data, 
for  example,  stator  and  rotor  fluid  turning  angles  and 
W2/W,  ratio,  are  easily  derived  from  the  basic  velocity 
triangle  data. 

For  purposes  of  comparison,  sets  of  velocity  tri- 
angles can  be  computed  at  the  same  streamline-plane 
locations,  but  using  the  unmodified  simple  radial 
equilibrium  axial  velocity. 

3.3   Computer  Implementation  of  the  Model 
3.3.1   General 

The  iterative  nature  of  the  detailed  calculations 
virtually  precludes  hand  calculations  of  the  streamlines 
as  a  practical  matter.   A  computer  program  to  accomplish 
the  calculations  has  been  written  in  FORTRAN  IV  language 
making  use  of  the  model  easy  and  inexpensive. 

The  goal  of  the  work  described  in  this  thesis  was 
to  develop  the  model  and  to  program  it,  so  that  it 
could  be  applied  in  future  design  studies.   Since  the 
model  itself  is  experimental,  the  program  is  not  merely 
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the  mechanization  of  standard,  accepted  practice.   For 
this  reason,  it  was  thought  advisable  to  include  here  a 
sufficiently  detailed  description  of  the  computer  program 
so  that  the  model  used  is  clear  and  so  that  changes  can 
be  made  easily  in  the  future  if  desired.   References  7 
and  8  were  useful  as  prograniming  language  guides. 

The  programming  philosophy  was  to  simplify  the 
understanding  and  use  of  the  program  and  to  provide 
flexibility.   The  program  was  segmented  into  subprograms, 
each  having  a  specific  function.   Identical  variable 
labels  have  the  same  meaning  in  all  segments  of  the 
program  in  order  to  simplify  understanding  it.   The  use 
of  COMJMON  reduces  the  storage  demands  of  the  program. 
Extensive  execution  information,  such  as  the  number  of 
loops  performed  and  intermediate  values  can  be  printed 
out  at  the  user's  option.   Diverging  iterative  loops  are 
terminated  by  internal  testing. 

The  following  paragraph  lists  the  significant  steps 
in  the  execution  of  the  program  to  modify  a  simple  radial 
equilibrium  design  in  accordance  with  the  model.   Sub- 
sequent paragraphs  will  discuss  each  step  in  detail. 
3.3.2   Procedural  Steps 

The  following  list  contains  a  list  of  all  significant 
steps  in  the  program.   Some  details,  such  as  initializa- 
tion of  variables  and  testing  for  loop  convergence,  are 
omitted  and  are  self-explanatory  or  can  be  conveniently 
discussed  separately.   Free  use  will  be  made  of  the 
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notation  employed  in  the  computer  program  itself.   The 
list  is  intended  to  serve  as  a  guide  for  the  discussion; 
full  appreciation  of  the  interrelationships  can  be  gained 
only  by  studying  the  program  flow  charts. 

(1)  Read  and  process  input  data 

(2)  At  the  first  plane  calculate  a  set  of 
streamlines  satisfying  continuity  using 
the  simple  radial  equilibrium  axial 
velocity  corrected  for  flare 

a)  the  first  streamline  follows  the  hub 

b)  guess  at  the  location  of  the  next 
streamline  and  iterate  until  continuity 
is  satisfied  between  tV';o  streamlines 

c)  find  locations  of  succeeding  streamlines 
in  turn  until  the  last  one  has  been 
calculated 

d)  compare  the  radial  location  of  the  last 
streamline  to  the  desired  tip  diameter 
at  the  plane  under  consideration 

e)  if  a  discrepancy  exists,  modify  the 
"average"  axial  velocity  at  that  plane 
and  repeat  from  a)  above 

(3)  Repeat  (2)  for  planes  2  and  3 

(4)  After  a  set  of  streamlines  is  calculated 
for  all  three  planes,  recalculate  them  all 
as  before  but  using  the  actuator  disk 
correction  on  the  axial  velocity  as  well 

48 


(5)  Test  to  see  if  the  streamline  positions 
for  plane  three  have  changed  much;  if  they 
have,  go  back  to  step  (4)  and  recalculate 
the  streamline  positions  for  all  three 
planes 

(6)  Calculate  and  print  velocity  triangle  data 
for  each  streamline  position  at  each  plane, 
using  the  modified  axial  velocity 

(7)  Repeat  (6),  but  using  the  simple  radial 
equilibrium  axial  velocity. 

3.3.3  Overall  Program 

The  overall  program  is  represented  by  the  flow 
chart  in  Figure  3-2,  on  the  next  page.   The  circled 
numbers  provide  points  of  reference  for  more  specific 
flow  charts  in  later  figures.   The  initialization, 
incrementing,  and  printing  of  LI,  L2 ,  and  L3  are  included 
in  the  flow  charts  so  that  the  user  can  readily  interpret 
certain  diagnostic  output  to  be  described  later. 

3.3.4  Reading  and  Processing  Input  Data 

The  first  step  in  the  program  is  to  read  a  set  of 
input  data  and  perform  some  preliminary  calculations. 
This  is  accomplished  by  the  main  program,  hlAlN. 

The. following  is  a  list  of  the  input,  determined 
from  a  simple  radial  equilibrium  design,  which  must  be 
submitted: 

(1)  working  fluid,  specific  gas  constant,  and 
specific  heat  ratio 
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(2)  turbine  mass  flow  rate  and  rpm 

(3)  total  temperatures  and  pressures  at  all 
three  planes 

(4)  hub  and  tip  radii  at  planes  one  and  three 

(5)  stator  and  rotor  widths 

(6)  coefficients  of  terms  in  the  simple  radial 
equilibrium  relation  for  tangential  velocity 
at  each  plane 

and       (7)  a  reference  velocity  and  radius  used  to 
find  the  simple  radial  equilibrium  axial 
velocity  at  each  plane. 

In  addition,  the  number  of  streamlines  desired  and 

certain  test  values  are  read  in. 

The  following  calculations  of  a  preliminary  nature 

are  performed: 

(1)  test  to  ensure  the  overall  flare  is  not 
zero  and  calculate  flare  angles 

(2)  calculate  the  basic  parameters  used  to 
apply  the  flare  correction  to  the  axial 
•velocity  (paragraph  3.3.9) 

(3)  calculate  dimensional  test  values  from 
input  (paragraph  3.3.13) 

(4)  calculate  blade  heights 

(5)  calculate  certain  other  often  used  combina- 
tions of  input  which  remain  constant  for 
the  data  set  under  consideration 
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Selected  parameters  are  printed  out  prior  to  further 

execution  of  the  program. 

3.3.5   Calculating  a  Streamline  Location 

The  heart  of  the  whole  program  is  the  calculation  of 
streamlines;  the  most  basic  step  in  the  strec'imline 
generation  is  the  calculation  of  the  "next"  streamline 
location/  RNEXT,  from  the  position  of  an  adjacent  stream- 
line, R(I,J)/  and  values  VXZ  and  VTZ,  representing  average 
values  of  the  axial  and  tangential  velocities  in  the 
channel  between  RNEXT  and  R(I,J).   This  calculation  is 
performed  by  subprogram  RNEXT (I),  whose  flow  chart  appears 
in  Figure  3-3. 

The  total  velocity,  V,  is  determined  from  the  axial 
and  tangential  components,  and  with  equation  (2.13) 
determines  the  static  to  total  temperature  ratio.   The 
temperature  ratio  and  the  equation  of  state  for  a  perfect 
gas,  equation  (2.10),  are  used  in  another  compressible 
flov;  relation,  equation  (2.17),  to  find  a  density  RHO. 
Then  the  continuity  equation  (2.35)  can  be  used  to  find 
a  flow  area,  AREA.   From 


RNEXT  =  /MEA   ^  R(i,j)2  (3.1) 


the  location  of  the  next  radially  outv;ard  streamline  can 
be  found. 
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Since  RNEXT  depends  upon  average  values  of  axial 
and  tangential  velocities,  which  are  in  turn  functions 
of  radius,  iteration  is  required.   This  iteration  is 
accomplished  by  the  subroutine  STREAM  shown  in  Figure  3-4. 

STREAM  utilizes  RNEXT  as  a  guess  for  the  next 
streamline  and  then  calculates  the  average  velocity 
values.   The  average  values  VXZ  and  VTZ  are  then  used 
to  calculate  a  nev;  RNEXT  and  so  on.   When  finally  the 
difference  between  any  RNEXT  and  the  previous  RNEXT  is 
small,  STREAM  returns  the  values  R(I,J+1)  =  RNEXT, 
VX(I,J+1),  and  VT(I,J+1)  to  the  main  program.   The  index 
L3  counts  the  number  of  times  RNEXT  is  calculated  during 
the  calculation  of  a  set  of  JJ  streamlines  at  any  plane. 

In  addition  to  finding  an  R(I,J+1),  STREAM  also 
applies  the  flare  and  actuator  disk  corrections  to  the 
unmodified  (simple  radial  equilibrium)  axial  velocity. 
The  corrections  are  accomplished  by  multiplying  the 
unmodified  axial  velocity,  calculated  in  subprogram  FVX(I) 
(paragraph  3.3.6),  by  the  corrections  for  flare  and 
actuator  disk  .effects  which  are  calculated  by  subprograms 
FLARE (I)  and  DISK (I),  respectively.   (See  paragraphs 
3.3.9  and  3.3.10)   Since  the  very  first  time  a  set  of 
streamlines  is  calculated  at  a  plane,  DISK  is  undefined, 
only  the  flare  correction  can  be  utilized  in  the  first 
(Ll=l)  loop. 
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3.3.6  Calculating  V   and  V^  (Simple  Radial  Equilibrium) 
In  a  simple  radial  equilibrium  preliminary  design, 

it  is  standard  procedure  to  pick  a  functional  form  for 
the  radial  variation  of  tangetial  velocity.   Then,  the 
radial  variation  of  the  axial  velocity  is  determined  by 
the  simple  radial  equilibrium  equation  (2.23).   In  the 
present  case,  the  tangential  velocity  can  have  the  form 
of  equation  (2.27)  where  the  coefficients  a_,  a^ ,  a  ,  and 
b  are  chosen  by  the  designer.   The  simple  radial  equili- 
brium equation  has  been  integrated  for  this  form  in 
Appendix  D,  resulting  in  equation  (2.30).   The  function 
subprograms  FVX(I)  and  FVT(I)  calculate  the  sim.ple  radial 
equilibrium  values  of  axial  velocity  and  tangential 
velocity,  respectively,  given  an  argument  RR  and  input 
parameters  A2(I),  A1(I),  A0(I),  B(I),  VX0(I),  and  R0(I). 

3.3.7  Generating  a  Feasible  Streamline  Set  at  a  Plane 
The  main  program  calculates  a  set  of  streamline 

locations  at  a  plane  starting  with  the  hub  and  using 
STREAiM  to  step  outward  to  the  tip.   Because  the  axial 
velocity  has  been  modified,  however,  the  average  axial 
velocity  is  different  from  the  simple  radial  equilibrium 
design  value,  and  the  last  streamline  does  not  fall  on 
the  fixed  tip  radius.   To  correct  this,  a  reference 
axial  velocity,  VXR(I),  used  by  FVX(I),  is  adjusted  and 
the  streamline  set  is  recalculated. 
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The  iteration  is  continued  until  the  last  streamline 
location  is  close  enough  to  the  desired  tip  radius.   A 
flow  chart  of  the  section  of  MAIN  that  generates  a  set  of 
streamlines  that  fit  the  geometry  is  shown  in  Figure  3-5. 
The  variables  S  and  SI  do  not  play  a  direct  role  in  the 
streamline  generation  and  will  be  explained  in  paragraph 
(3.3.11).   The  adjustment  of  the  reference  axial  velocity 
is  explained  in  the  next  paragraph.   The  index  L2  counts 
the  number  of  iterations  made  at  a  plane  to  converge  to 
the  proper  R(I,JJ).   At  the  users  option  (by  specifying 
KK=1  in  the  input  data)  L2  and  other  parameters  may  be 
printed  out  after  each  trial  set  of  streamlines. 
3.3.8   Adjusting  VXR(I) 

Initially,  VXR(I)  is  set  to  VX0(I),  the  value 
deterroined  from  the  simple  radial  equilibrium  design. 
Whenever  the  location  of  the  last  streamline,  R(I,JJ), 
misses  the  tip  radius,  RT(I),  VXR(I)  is  adjusted  to 
increase  or  decrease  the  average  axial  velocity  before 
recalculating  a  new  set  of  streamline  locations,  in  order 
to  make  R(I,JJ)  converge  to  RT(I). 

To  find  an  expression  for  a  correction  to  the 
reference  axial  velocity  in  terms  of  the  error, 
DR  =  R(I,JJ)  -  RT(I),  an  expression  for  the  flow 


ni  =  V  -A  (3.2) 

P     X 
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can  be  differentiated,  giving 


d{|)  =  AdV^  +  V^dA.  (3.3) 


But  we  desire  d(m/p)  =  0  and  A  =  Tr(RT(I)^  -  R(I,1)^) 
whence 


2V  RT ( I ) • DR 
RT(I) 2-R(I,l) ^ 


From  experience,  to  make  the  convergence  more  rapid,  the 
following  term  was  added  to  the  VXR(I)  by  MAIN: 


^w   =  2  VX(I,JJ)  DR 

^^x    RT(I)  -  R(I,1)  ^    '' 


3.3.9   Correcting  Axial  Velocity  for  Flare 

The  simple  radial  equilibrium  axial  velocity 
calculated  by  subprogram  FVX(I)  is  corrected  by  multiplying 
it  by  a  non-dimensional  factor  computed  by  subprogram 
FLARE(I).   The  basic  relation  involved  is  equation  (2.33), 
rewritten  as 

FLARE  =  cos^cf)  =  ,  (3.6) 

tan  ^4)  +  1 

where 


tan  <P   =  ^^  p  ^^^.  (3.7) 
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RR  is  the  radial  location  of  the  point  of  interest,  D  is 
a  parameter  depending  upon  the  plane  under  consideration, 
and  REF  is  a  geometric  parameter.   For  additional  details 
see  Appendix  F.   FLARE  is  always  positive  and  less  than 
or  equal  to  one. 

3.3.10  Correcting  Axial  Velocity  for  Actuator  Disk 
Effects 

The  actuator  disk  correction  is  applied  after  the 
flare  correction  and  is  also  a  multiplicative  non-dimen- 
sional factor,  calculated  by  subprogram  DISK (I).   The 
relations  used  in  DISK(I)  are  forms  of  equations  (E.l), 
(E.5),  and  (E.2)  for  planes  one,  two,  and  three, 
respectively.   DISK  is  always  positive  and  may  be 
greater  than  one. 

3.3.11  Testing  for  Overall  Streamline  Convergence 

The  corrected  axial  velocity  at  a  streamline  depends 
upon  its  location  and  the  value  of  axial  velocity  on  the 
same  streamline  at  planes  up  and  downstream.   Therefore, 
in  general  it  is  necessary  to  recalculate  the  streamlines 
at  each  plane  several  times,  the  streamlines  moving  less 
with  each  iteration.   The  streamline  positions  are 
calculated  in  order  at  planes  one,  tv/o,  and  three  at 
least  twice.   Each  time  the  streamlines  for  plane  three 
are  recalculated,  the  following  sum  is  accumulated: 

JJ 
S  =   2   |R(3,J)^j^j^g  iteration  ~  ^^^'^^last  iteration 

(3.8) 
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Refer  to  Figures  3-2  and  3-5.   If  S  is  small  enough,  the 
streamline  positions  have  not  changed  much  from  the 
previous  iteration.   This  means  that  the  set  of  corrected 
axial  velocities,  VX(I,J),  and  streamline  positions, 
R(I,J),  have  converged  to  mutually  consistant  values, 
and  iteration  can  be  terminated. 
3.3.12   Velocity  Triangle  Calculations 

Once  the  streamlines  have  been  located  at  each 
principal  plane,  a  velocity  triangle  can  be  constructed 
at  each  point.   PI,  P2 ,  and  P3  in  Figure  3-1  are  the 
three  points  for  one  particular  streamline.   JJ  triangles 
are  constructed  at  each  plane. 

The  fundamental  input  variables  to  the  velocity 
triangle  calculations  are  the  streamline  location  R(I,J), 
VT(I,J),  and  VX(I,J),  which  were  calculated  concurrently 
in  the  iteration  process.   The  following  derived 
quantities  were  calculated  by  MAIN  at  each  plane: 

U  =  2TT  •  RPM  •R(I,J)/6  0  (3.9) 


V  =  /VT(I,J) 2  +  VX(I,J) 2  (3.10) 


W  =  /VX(I,J) 2  +[U  -  VT(I,J) ] 2  (3.11) 


a   =  arctanC^"-^^-^"^^  )  (3.12) 

a  -  arc-can  ^^^^-^^jj  j 


VT(I.J)  -  U)  (3.13) 


^  =  arctan(   vx(I.J)   ^ 
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M  =  from  equation  (2.16) 


M   T  =  M  •  5^  (3.14) 

rel       V  \  J .  -LI ; 


In  addition,  the  follov;ing  quantities  were  calculated  for 
each  streamline: 


VI  =  V (plane  2) 

V2    V(plane  1)  lJ.ii>; 


W2  ^  W (plane  3) 

va   W  (plane  2)  U.ib; 


stator  turning  angle  =  a  (plane  2)  -  a (plane  1)  (3.17) 
rotor  turning  angle   =  3  (plane  3)  -  3  (plane  2)  (3.18) 

To  compare  the  above  results  to  those  attainable  by 
simple  radial  equilibrium,  all  of  the  above  velocity 
triangle  calculations  are  repeated  at  the  same  locations, 
but  with  the  unmodified  simple  radial  equilibrium  axial 
velocity. 

Figure  3-6  is  a  flow  chart  describing  the  calculation 
and  printing  of  the  above  results. 
3.3.13   Iteration  Test  Values 

For  flexibility,  the  user  supplies  the  test  values 
used  to  terminate  the  three  iteration  loops.  The  three 
loops  are  (1)  overall  streamline  convergence. 
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(2)  convergence  of  the  J J ' th  streamline  at  a  plane  to 
RT(I),  and  (3)  the  location  of  the  next  streamline 
R(I,J  +  1).   These  three  loops  are  counted  and  labeled 
by  the  variables  LI,  L2,  and  L3,  respectively.   Test 
values  TESTl,  TEST2,  and  TESTS  are  calculated  from  input 
parameters  Tl,  T2,  and  T3  and  the  following  relations: 

TESTl  =  Tl  (RT(2)  -  R(2,l))  (3.19) 

TEST2  =  T2  (RT(2)  -  R(2,l))  (3.20) 

TEST3  =  T3  (RT(2)  -  R(2,1))/JM1  (3.21) 

An  iteration  loop  is  terminated  normally  when  the  error 
term  becomes  less  than  the  appropriate  test  value. 

3.4  Convergence  Checks 

If  for  some  reason  the  program  ever  begins  to 
calculate  diverging  values  of  the  error  term  in  any  of 
the  three  iteration  loops,  the  program  will  automatically 
terminate  the  calculations  on  the  current  data  set.   This 
is  accomplished,  as  depicted  in  Figure  3-7,  by  exiting  if 
the  error  term  has  increased  at  least  twice  in  that  loop. 

3.5  Program  Verification 

In  this  thesis  a  model  to  improve  simple  radial 
equilibrium  design  by  applying  flare  and  actuator  disk 
corrections  to  axial  velocity  is  described  and  a 
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computer  program  to  implement  the  model  is  presented. 
The  computer  program  has  been  verified  to  the  extent  that 
there  is  confidence  that  it  accurately  reflects  the  model. 
Full  evaluation  of  the  model,  including  comparison  to 
more  complicated,  exact  design  methods,  must  be  left  for 
future  efforts. 

Although  the  iterative  nature  of  the  streamline 
calculations  makes  hand  calculation  of  the  whole  procedure 
very  burdensome,  it  is  easy  to  verify  a  solution  by  hand. 
The  outputs  of  subprograms  FLARE(I),  DISK(I),  FVX(I), 
and  FVT(I)  were  all  verified  by  hand,  given  a  computer 
calculated  streamline  position  R(I,J  +  1).   With  these 
values  and  R(I,J),  RNEXT(I)  and  STREAM  were  verified  by 
hand  calculating  R(I,J  +  1).   The  "computer  calculations 
of  flare  geometry,  velocity  triangle  related  output,  and 
other  algebraic  relationships,  all  agreed  with  hand 
calculation. 

Since  each  input  data  set  must  be  a  complete  simple 
radial  equilibrium  preliminary  design,  the  number  of 
independent  geometries,  working  fluids,  tangential 
velocity  distributions,  etc.,  tested  was  not  very  large. 
Nevertheless,  it  is  felt  the  verification  of  the  program 
was  sufficient  to  ensure  accuracy  for  feasible  designs. 
Results  are  shown  for  three  sample  cases  in  Chapter  4. 
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4 .   RESULTS 

4.1  General 

This  chapter  presents  the  results  of  applying  the 
computer  program  to  three  sample  simple  radial  equilibrium 
designs.   The  computer  output  listing  is  included  and 
plots  of  representative  variables  are  provided  for  easy 
comparison  between  the  simple  radial  equilibrium  input 
design  and  the  modified  design. 

4.2  Free  Vortex  Design 

Table  4-1  shows  the  characteristics  of  an  axial 
inflov;  and  outflow  free  vortex  stage  designed  as  the  last 

stage  of  a  double  flow  turbine  in  a  geothermal  power 

9 
application.    A  sketch  of  the  geometry  appears  in 

Figure  4-2.   Notice  the  design  has  a  constant  hub  radius, 

but  has  large  flare  in  the  shroud.   The  working  fluid  is 

Freon  21.   Eleven  streamlines  are  called  for. 

The  computer  printed  results  are  listed  in  Tables 

4-3/  4-4,  and  4-5.   It  is  seen  that  the  large  changes  in 

axial  velocity  (about  25%  at  the  hub)  can  result  in  a  six 

degree  change  in  the  relative  inflow  angle  to  the  rotor. 

Less  pronounced,  but  readily  noticed  changes  are  shown 

in  other  flow  angles,  velocity  ratios,  and  Mach  numbers. 
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In  Figure  4-6,  the  axial  velocity,  stator  outlet 
angle,  rotor  fluid  turning  angle,  and  relative  velocity 
ratio  are  plotted  for  both  modified  and  original  simple 
radial  equilibrium  designs.   It  may  be  significant  to 
note  that  the  changes  at  the  hub  (where  (p   =    0°)    may  be 
considered  to  be  "indirect"  in  that  they  result  from  the 
increased  average  axial  velocity  required  to  obtain  the 
same  mass  flow. 

4.3   Constant  Reaction  Design 

The  input  parameters  for  another  example  having  the 
same  annular  geometry  as  above  but  with  a  different 
assumed  tangential  velocity  distribution  are  shown  in 
Table  4-7.   Axial  inflow  to  the  stage  is  maintained, 
but  the  constant  reaction  type  simple  equilibrium 
relations  (see  Appendix  D)  are  used  at  the  rotor  leading 
and  trailing  edges. 

The  computer  printed  results  are  shown  in  Tables 
4-8,  4-9,  and  4-10,  and  plots  of  axial  velocity  and  rotor 
relative  flow  angles  appear  in  Figure  4-11.   From  Figure 
4-11,  the  modification  to  axial  velocity  doesn't  seem  as 
significant  as  in  the  first  example;  the  influence  on  the 
rotor  flov;  angles  is  still  significant,  however 
(A6  =  6° "at  the  hub) . 
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4.4   Free  Vortex,  High  Mach  Number  Design 

To  demonstrate  the  adequacy  of  the  computer  program 
to  calculate  in  regimes  where  M  >  1.0,  the  preliminary 
design  summarized  by  Table  4-12  was  run.   In  this  case, 
the  flare  is  symmetrically  distributed  between  hub  and 
shroud,  as  is  shown  in  Figure  4'-13. 

The  detailed  results  are  listed  in  Tables  4-14, 
4-15,  and  4-16.   The  axial  velocity  and  Mach  number  at 
the  rotor  inlet  (plane  2)  are  plotted  in  Figure  4-17 
for  comparison  with  the  simple  radial  equilibrium  values 
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tv  if.  >c  in  fv,  vc  a. 

tv.  s  >o  s.  M  a  ^ 


ifi  rv.  ^ 
>£■  <c  a, 


cv,  ^  ^  I    I   ^  r-  a  rv 

I    I    I    I 


•M     t.   "^    r    i    "^1  "}    Si    rVj   J-;    J- 

r"f>.(rr^ojCj.(S..-i.»a  ^£ 
ajnjif)r«0"^r,  if.  r^Sy 

N£,.c-c-o>c<;!vr^r.rvx 
•c  -c  -c  <.  <,  •c  <  <c<c>f<o 
t     I    I    I     I    I    I     I    I    I     I 


u.  >C 

•-  ^       ^J 

c 

o  X 

r  ■ 

z  — 


s.  •=;  s  -s,  •"»  s   s  -s  s  ■5.  ^ 

•iS5.t>.I«6.S.C»S5.&. 


i.  t.    ^i    C;« 


^i.   i> 


x-cif  iTxx  rvt-cr  a.  a 

\r.  s  <^  t^  <  r.  a   -t  ■*  ^  a 

r~   -r  X  <.  ri  J^  <.  ^  3  ii    a. 

X  r  -o  «  o>  r^j  (^  -c  r-v  .rH  rr 

r^  r>  K  J-  6.  r~-  tr.   ^  (T  ^^■  -c 

x^ttrir— vcaxj^^x 

fv  (»)  (^1  J-  J".  XI  -c  >o  r^  X  X 

(^  r-'  r"  (^  (^  rt  r-'  rr  r^  (^  c 


■5.    'Sti&^S.t.    S.5.S.t. 

S(S.S.SS<S    S.S.SSEy 
•C     <   <^     -J^    -C    <    ^    ^   <     -Ji   •£ 


•-<  .£)  ri  n  0^  — <  n  X  i  o>  r^ 
jiLDf^'-^Txojrucv'^X) 
X)<r  r^<.rfOX  .o.s  ojror^ 
x(r«j-,ti(rnx   nj%.c 

t\.  (c  r,  r.  t'i  r.  j  j^   j    x.  vt 


fofVJ^O"^£)x^^r^  a^  >c  >c 
&u".  m>OvDxnixixs.ir 
-»  x   'J  r".  i"^.  a.  r^  -»    J"i  <*  if. 

"fn<rr^.f>NrxwNC*>o 
cr  X'  rv  >o  -£  Xj  Xj  .*   j-  j-  -r 


<05>.S.*5iK(V.X    ■*—•<. 

ci  ■*   K&.>£ixr^x.  i^j^nj 

L'"/  X)  r^  Oi  ^  •*  r^  s  m  X)  X 
ca    r^.f'^cxi.t-r'^rv'— 1 


■^   x~  /i  'S.    Ti  %  r-.   ,C    X.  "i  rr 

<  .»   -i.  r^  •-•  >£.  a  X)  X)  rv  t>.  Oj 
I        X.  x<f^c».mrv(\j(sv(vicx,f^ 
a.         ........... 

_.       xj^j^NXNCa.moj'so" 

<  CV,CVj(VlfVj--l»-»*^,-i^rHi 

I     I    I    i     I    I    I     I     I    I 


X  o  n.  cs  re  f^  -t  •-•  ^c  X  r^ 
s.   (T'rHOjxtrics   yDr~.(\j-r 

0  r".  ryu,^o.-»  f^x.X'i/ 
S'  -»  X'  X)  X  ^o  *  j^  nj  rv  m 
frro.*NCF.»"SMrv^ 

3^  a  (s  s  •^  rv,  oj  D  n  J-  X. 


s  -c  .-  r».  ^  ci  X.  f^  .H  X  oj 
vc  a  X  •-  X  rv  a  e.  ai  -»  >c 
c.  of^-r-ci^.-.  •-.  Ci  c  rv. 
•-^xi^.rx.  ■£r~-X)crxr) 
r^a   t>r^i^<\jxxcvjs<r. 

■£i  ■»  .^mfvjfv— 1  — <^^^ 


{\J^cx)^c\J(^(^J(^J— "CTivC 
s.x.»x-»xr-.  ^'-x  — 
rv  -r  r^  r  rvi  r-  X  a  ir  rv  x 
xix  x/nfnxisi~-  <Cno>c 

•-•  r.  '^  ^  .»    ^  r.  r^  Cs.  p^  A. 


fv  Xi  X  <vj  -tj  f^  D  li,  r^  xj  X 
f^NCX.  mat^xma/rvo- 
ojx,  Xja-sfvax  a— it, 
ro^CsC«-ijN£>ru(r^.c(r 


X, 


^t^CX)XlJ--»     -»-DO0J 
llllllllll 


> 


&e>.SC>.!>.S.G.CSSS(S. 


X  tn  nj  (V  X(  -»  m  cr  x  rv  cv 
ccvr-o(s  eaa  -je^a 
o  ro^Cif^<^J^c^fnx)X) 
cv.j-o--t(s  vOX.-'-ccrtr 
(\ij\(^j^x.  -jojx   t.   x»- 

cixi6.xisx)s.»  a-  t\j  <i 
ar-r->c<xxj^rr-C- 


i   r^m.*   .»-«ix   x(n^ 
(T    .*    ax-*'f^(T    a    cu    a 

(Txiittxfv.— ifv    f^xirv 
xr.  r-tt    aox^c^Oi.<c 

oos.--«cvjn-»-X5vOivaiic 


cuocfvj^ans  4^*  «cr 
0"   f^>  •£  IS   n  X  ro  <-•  cu  XI  fT 

a.   a  a   &.«-••-■   fXjfor^;   *X) 
fViOjojr'ironaimmDD 


futvs>cx^tnr5crx>>* 
cr  -i    h-o^Sforv(»)Ojxn 

X  CTi  c  ti  -^  pj  D  ^    X.  >f .  r- 
cvj  cvj  n  ri  n  n  n  m  m  (^  fn 


<vi  r^  K  <o  c^  tn  (v;  «■  (r.  fVi  xi 
.»x>Jl^ox)•-■n^\Jxs.^o 
cr  r-cvjxr^Ein  xxir^nj 

j^xi^ONO-ct^r^r^xiiri^ 

xo^e;--'t\jr;  -^■X)-olvtJ^ 


c    <r  tf . .  r^  c  <v   a   x   o  .  f  »». 
e>  r:  X    <c  o  .»  n  <r  c^,  X  X. 

&.0.-»      .<;x!a.(\.X.    r^O.    fVj 

[7    C(r    ocf.    f^e,cs»- 
<a  s  IS  s:  t»  -•  "rt  .rt  .rt  .1  .rt 


K6.    .»-X;Xrv-»XC\jXO- 

r-*is>oajocx>f>J--i'S"-< 
6»  f,  .tx.-<nxo  (vx.  X 
OCT   o~<r^6.sis^--«'^ 

SIS»Sl^^^.-«'-<.-ir< 


s  D  ui  %<  s   X  tr    ^  xi  r   * 

SJ    OC-fCVjI^X     <S.0J.€J- 
rvw.»rvs    (VI  X   (rmxX) 

ao   (rfisb.w.^^(v<v 


•-■fVin.»X)-orvX(r>S"-< 


•-4(\jnj^x)-orvoo<r><a"-« 


•-<  fur>.*^x)-otvoc(n'Si»-< 
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RrSULTS    ra    SRp    MrolPIf^D    r\y    riiPE,     DIs< 


.1 

Vl/VP 

W?/wi 

DALP'^A 

I'tF  I  A 

1 

2  .  v*.  t  ■■  ?  1 

2  .  (7. 1  6  7  -3 

bl  .32275 

-h^.t,7•^^'' 

p 

2.;i?lr>i 

2M1371 

6  1  .  1  6  U- 1 

'Ih.iZZ?-) 

"» 

a.iM  H6 

2.2?ff^s 

61  •2iJ2?h 

-73.'»7>^-^t 

« 

?.??U2 

2.33?'.*, 

61  .6'5'='PX 

-7/.  16V  6'. 

c 

?#  cri.^i,'. 

2.'.71t7 

62.26S75 

-66. /93h7 

A 

?.C8=S>5 

2.f  157JJ 

63.-/^51? 

-k3.^?'<£'' 

7 

2.1'.la6 

2.768F5' 

6't..'  X5>6 

-e>9.23f  22 

R 

2.  21 '-'.9 

2.92?f-? 

fe5.2'»S6'^ 

-S-..  73bl'. 

q 

?.  l^f'^RS 

3.p6^r? 

66.551f-7 

-'♦9.63699 

it- 

?.'.27t7 

3.1777'; 

6  7.97316 

-43. ^*5/6 

It 

2.57U0 

3.2417? 

69.<,qiQ3 

-37.79771 

repilts   p-or   q^F,    iiN^ionif  ifo 


.1 

Vl/Vp 

w:>/i-l 

'''ALPHA 

■^Hi   TA 

i 

P.'.f"  >72 

2.''?2?6 

62.5;5Sf6 

-r6.'5l^"77 

•? 

?.35-'3'» 

2.  t89>.« 

6?.Q95^9 

-<<2.  378S<. 

•» 

p.??  177 

2.325'.S 

63. l???** 

-7«.''.°bb7 

<i 

?.?5''f6 

2.i,57(.1, 

63.'.i?,l?6 

-73.6281'. 

c 

?.?l?'Ac^f 

?.5836'' 

63.669Q? 

-6B.95/Kb 

t. 

?.  161'.9 

2.7P12R 

63.98«li? 

-6'.  .''.'»S75 

7 

2.11751 

2.SP71', 

ft**  .  362=:? 

-bfi.^t'-Sbb 

SI 

?.('7'.13 

2.R974  1 

6'».8^2>:  1 

-53.o2'.63 

Q 

2. ^3  "°2 

2.o6?^23 

65.32106 

-'.7.92t  93 

1-7 

1.987S? 

3.t?15?7 

65. 93703 

--.p.  11613 

ll 

1  .9',Tj,7 

3.J>37T? 

66.6757C 

-36.  13>'93 
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f: 


-76 


200 


X 


lliO 
(ft/sec) 


I'fodified  SHE 
•SHE 


V^ 


FiRui-^  h-11 
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RUK    NO 


SIMPLE    ^ADlAl      EOuiLlr-^IUf      ^^Kl.IM    Tl!R■^I^:E     OESlCiM 
WITH    r  L  A  P  r     A  \i  D     A  C  T 1 1 A  T  f  P     0  I  S  <     C  n  c?  R  c-  C  T  I  P  \  S 


G  A  M  "  A  =  1  •  1  3 P" ('"  H q  =  .^^P.^i'^- .'- ^ 

PGA?=  It"!  t  Pf'^vp'  HR=  0.361=^M 

RPm=  TQi;^.C?7f'  TrSTl=  B.C3!>^3?lf: 

FLOw=  ?36*6?7?  TESTa=  0.02lf>-^5 

REF=  f  .qv-^??  ?•«'*  TEST3=  7j  ,  ?'''l' \bl 

■^E=  r.a'tSQPS  TEST^=  0.;v7.t'r^e 

•^  R  =  r  t  0  8  3  ?  '7  y  w  R  =  /^ .  Z  ?  3  ^/  'J  P 

PMTlP=  26.ri?a  PHIHU3=  -26»"i:i 

JJ=  11 

^'LAr-E  1  2  3 

Ttki)  ?«8»0000  ^,88..'"vl0i-  bB'-'C*?''') 

P.?II)  60»0P)0C  6J'«v'3^('  33»6e^5B 

PT(I)  1  ♦026999  l.C675:^C  1  .  i  .-^FZ  J^'!5 

R(l,i)  p.y^^yc^^v"*  p.7465?e  (?. .  7  ^e'/.^K 

R:-)!!)  C^  .9''a7C^0(?  ?*Sv)77y7.v  2*?''.7?/jd 

VX0(I)  133»OC»0^  ?i^.vi/5?-v?  13'5.f''0.-;0 

B(I)  ['Uv^?^?!^?  ^('^^•l?99r  e.*i?0V'0 

A0(I)  C^j.00p"0r:'  0»'/10??C^  P.  ^00^50 

Ai(:)  0.00000  0t  0000(7  0.000^0 

Ap(I)  0.00000  0. '10000  0./0000 
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stator 


-^ 


tit) 


hub 


■< rotor 


^h,t=  26.0^ 


r^  =  0,907' 


axis 
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^   i  .*v    f>  .s.  s   -s  .-»   i.    i,   -. 


^  .-1  r,  N  ^  r^  c\.  c^  f^  ^  i 
•-'  ir-"SXErvf^r^xs.-»- 
-^       o  <x.  r^  ^  <:  .r  s  -r.   ^  -t  -t 

i.    r.    ^   C:.  ^    ^  ::^   '.    ^  t.  ^. 


J".  X  -»  r-  r^  r^  (V  u"'  rv  r.  3^ 
TO  a.'rvvc.  -^soc  in.-'r^  * 
c^r»xaccx?~.  inovb. 

1      ^.    -u    »i     ».    -.    •-.      •>     -i     V.     rt 


T  <<  e\v  <\}  a  f  o.  w  a  ', 

X  :»;  rr  a  ir.  t.  ^  rj  ri  r^  -r 

-"  J  •-  •^  <  <  ^  r^  •-  ii  c 

<vi  <M  (Vj  .->,•  fvj  (\-  fv  fvi  a  (V  (v 


■i    N 


«i     5.     S 


■xn;:fj~^(r'^r-    x-»v£ 
cr  ;>.   tv  m  ;;.   ^  t.   ^  ;_  X  ^ 


r^    f^n.'j^    .-ttrn  ^cc^    (T- 

X    X   -»  c    (\r  t.  J  s    vC  s.    r^ 

c     i    J    J    n    ■<  iv  f.  ; .  <    r 

^t^rrr-.  fp  trar.f* 

'VI    ru  f\'  r^    pi^  rn  m  fvi  (\)  f\'   (Vi 


ir  •-"  — •   X  f".  j^    cr  r^  ^   r^  (\. 

c>   -*■  0"    N.  (T.  s    i^  ^  r)  Lf'  *-" 

ecr.  xxxo-crcinjr^j- 

I    I  I     I    I    I     I    I    I     I    I 


x^   no-p^-Lxrvxixrv 
p~  tr   u~   r^  J^.  -c  r.  -r    ti  .-■  if. 

^       fT)  y-   J-   (T^  X  r:  J'.  tr«   -^  ,c  s, 


iT  a    rv    ■  r^  rv  "X'  iTi  '-  ■'^  tv. 

f^  oc  sc  Oj  X  -c  ir.  f^  •-'  ^c  » 

>t  i  IT.  X  N  rv  r-.  fj  J-  c.  c^ 

5-  ir-»i^-»fv.minoccvjir: 

X  t^i^r^  xcr.  t.  ^(\,-»-  ji 

nC'  'ijNC^f:  ^£<^^  r^rvr^  r^ 

I  I     I     I  I     I     I  I     I     I  I 


&(s>.sj.%ej6.cs.ss.<sfs. 


S.    (S   C- 


C    fi    S    CS    Cs 


•-"       rvjifistvjriS'^'Htvj'-'.fi 

^  oj  vo  :r^  'H  -5  -s  X  X!  •-■  cri 

<        (ri;i.crxrv^a>ci-£itvjD 

I        .ocr\x-rir)S.aosim(r.c 


J^   f^    .i)    ^o   -c 


X  f^    rv  r^   X 

-c  -c   -c  -o  -c 


S    S   S  S   ri 


?S.  -2    !i  Pi  -^    S; 

s&.ri.c>s.si(s:siti 


n  ^  .£  •£  ^  o  vc  X  -»  rv  -» 

■r  C\  r;  r^  cc  o  s  (v.'  oo  x  x 

-i,  -r  :^  ri  ^  r^  -»  c\.  I.'/  (\j  ai 

cr^x  ^x-?  ^r^rv.'  •-•;») 


0J&  r)S--'(\.(\aO.G.>c>-' 

s  fvj  fxj  \C  C'  -^  m  0\  cu  &  cr» 

ii.    •-■    CV.    iT.  J     -L    ;i.    fVi    r/  rj   iT; 

f^xxxrrojxr^sxoc 
t-^crojO'Sirirx  -rr'^ 

,-ioro  ^(Tf^j^n  (vjsx 
-r   r-p-r50jcvj(\jOjfvjaj.-' 


r^  (\.  -»  cr  rv. 

^  rf  ^  r\j  -J- 

*  c.  iT,  m  -c 

(\j  <r  •-'  r^  rv 

»£>  >c  in  m  E 


r^  r^  X  -*  'O  ^ 
rvj  m  X  s  nj  cr 


X  ^  Kj  X  cr —-  ^  X,  (^.  c  ;r 

■f<C>XOu.— >t.t^fv.o  c^ 
%c  fu  -r  r^  (\,  (\j  iT)  oj  r^  -J.  iT. 
IT,  r>-*^tviri^iri'^fv.m 

mr^if)s.tviOj<sir)snir) 
^  n.'  m  -»  J-  J   -f  n  n  ,-v.'  — 1 


ir  a.  -t>tax5%j^t.c\.rL 
.£S  ^cm^■.-•tVJr^  -j-rrin 
^  4-   00  X  'S.   -»   »  <\J  r-  (\j  f~- 

(^>Cri^^XyDlflfT)(\jS 

rvso  vt-o-oioifiir)  inin'O 


r^   ir   r^  fn  c  r^  ^  IT.  -»  -»    vO 

xxr-t-t.P'/^a>cv-c>t. 

<r.   cr.  >c-T--iO^'^rniy(VjX 


'i    '-    '^     S   s    ^     •     "     "v    -^   -s 


t.  ^ 


^..   v:. 


lS-»a>CXXC:.ir.     .Iff,   Ci  Uj 

m>-•OJ>C<lC0OXM^C•-•^C  z 

--1        IT  cv'  c    •"  u    f^'  -»  3".  'v  %r  --I        _:       ."v 

ivai.ijr<t\j<r>crvs&i>^       a. 
I-        ...........  I- 

>        n£  --I  o   IT.  J.  r-  Oo  c    -fi  -I  a,  > 

r^rocrir^iDrocvjii   o^oor^ 
>c>oirioinirinir  j^^* 


S     (SSSSSISISSISS 

;    ^    «    •    -    "    ^    -.  X   s   -; 


(\j        if<sooi.«.-i>rv(\j<r.  s 

<f  (V  -1-  (^  fv,  a   IT  tvj  r^   s  in 

>        iff»;^jr^ir—  a—  fv:f 


If 

rv 

ir. 

s 

<\j 

Ol 

Q 

in 

« 

"5 

f. 

in 

- 

0. 

-» 

-t 

-» 

-t 

r 

fr 

a 

- 

K 

in 

.-* 

f^  IT) 

r» 

CT 

S 

o 

X 

-» 

•0 

tn 

-» 

m 

w-t 

^ 

<c 

r^ 

■X 

»-• 

• 

a 

a, 

-*■ 

00 

• 

(Vj-»ir)in-»--»  Dfu--«— <^ 
jvjm^   in»^r^  x<riCi..-'Oj 


(VX^vO-rr^-fr^ajnr^ 

f\.-t  inir.  r^-»  r-i-~  (^v-HO^ 
irocs  -»->ci^(vxcvirr~ 
(X   (X  ir   r^  •-•  >£  fs  X    X  m  X 

ifi-ca-cir.  i^nin  >c.£>-* 

^tKj(\J    (V/CVJfMfVJ<\jf>J--i^ 


r^r^in(rr^Gxx  \C--«(Vi 
If  1  Ci.  r5  X)  r^  -»  X.  in  oj  n  <r> 
-c  tv  cr  in  ir  If .  r^  if.  s  n  iT 

inrvj^  cnoj^in-orvr^f^ 
soj-*inr-x   0"is.-'(\jfn 


r^    P)D-*^-«xr^ini>.n-c 
6.     t    c-  if    »-■  6.    vf    6    f~-  J    •-■ 


-»  *rv(\jinnsom-»'-'<vj 

s  e,.rHip,os(ri-o<^j"-<(\j 

(T-  ^s.r»ncrin«noo   (^ 

X  ^p)-»   >or^(rscvjmin 


S    .-E>'<iM-J    -»    &.a    Clf. 
s^cr>.oc<VJCT■f^Jlr<c^(J^^~• 

X--J   <£cr'-inifx5>c\j 


c  ir  f  -»  J  If/  f\.  ro  Ol  >&■  <v 
If  r")  ^  .^  tr  — I  vc  o-  j^  .£■  r». 
'OS  fvj   xir>crt  ^ooJxr^oc 

-JCv-tl^CP,    O0-.    .--I^C 

r^x«xcraaa(s<ss 


s  i^^r^^-c-H  !sKcr  oc 
>o  (T'NCxxtrtxjocfocr^o 
^  ins-f  xf^.-^;an^cci. 
r^-rvxa    xcr(r<r.ss— • 

S     S.S5ifiBlciSl.-«^.-< 


^-ff\Jf^^^n^Of^ooJ^^*^•-< 


~j       .-i(\jm*in>orvoo(r\G(. 


~i      .-tcvim-j-invor^oocri'^"-* 
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t^  .>  -c  *  r^  rv  i  4,  s.  X  O! 

(T.  ff>  (T|  J^  >£,  f^  ^  x>  •^  (V  IT 

C.  »v  X  N.  ,C  ^  fV  -.<  JP  rv  X 

IV  rv  r^  r^  ,x  JL  -x-  <  X  X  r 

:-  "^  ■"-  »  i  Tv  .^  :l  ^  .S.  s 


•t    ^    H.     i    ^    X     -t    C^,     -L     C    -t 

c  ix  IV  ,r  .c  '/■  .n  iT    J   J  ^ 


^  :s    s 


X      ^  »-.  ^    i   X  J,    ^   r\.    ■.     X 

JK    nj^crc^<:c^fvxx»- 


^:    r    s 


fv    (\    ^    Tvj   fv;    tXj  (M   o.   !\j    (V    (V 

>;    >i   vf    \L   M   vi  »f   -c    r   <    c 

rvrv»vrvrvrv  fvfv'vrvrv 

iV    vj  .X  .X  (\.  (V,  (V  a  ru  Cv  r.. 


N 
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5.   CONCLUSIONS 

5.1   General 

A  model  for  improving  a  simple  radial  equilibrium 
design  by  modifying  the  axial  velocity  has  been  discussed. 
The  computer  program  presented  here  provides  a  means  to 
quickly  and  inexpensively  incorporate  the  model  into 
preliminary  design  procedure. 

The  improved  preliminary  design  model  is  basically 
a  perturbation  of  the  simple  radial  equilibrium  design 
procedure,  and  comparisons  of  the  two  procedures  can  be 
made  easily  from  the  computer  program  output.   It  would 
be  useful  to  compare  the  improved  model  results  to  designs 
produced  by  more  rigorous  and  complicated  analytical 
procedures  as  well,  to  see  where  the  balance  of  complexity 
versus  accuracy  lies.   Since  the  program  described  here 
is  only  a  "design  tool"  and  was  not  derived  analytically, 
comparisons  to  more  exact  procedures  are  really  needed  to 
give  full  respectability  to  the  results. 

The  changes  made  as  a  result  of  applying  the  model 
to  a  simple  radial  equilibrium  preliminary  design  are  of 
a  magnitude  which  suggests  that  they  may  have  measurable 
effects  upon  turbine  performance.   Further  study  is 
needed  to  quantify  these  effects  in  terms  of  efficiencies. 
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5.2  Improving  the  Model 

Aside  from  possible  improvements  in  the  generation 
of  the  correction  to  the  simple  radial  equilibrium  axial 
velocity,  other  changes  in  the  program  to  extend  its 
usefulness  may  be  desirable.   An  example  is  extending  the 
range  of  the  program  to  include  more  than  one  stage, 
with  two  conical  sections  per  stage. 

5.3  Applications 

The  most  direct  application  for  the  program  is  to 
refine  the  blade  shapes  of  an  axial  turbine  in  preliminary 
design.   However,  less  specific  applications  to  research 
type  problems  may  be  fruitful.   Used  in  conjunction  with 
a  sufficiently  general  simple  radial  equilibrium 
preliminary  axial  turbine  design  program,  the  program 
presented  in  this  thesis  could  be  used  to  generate  a 
series  of  modified  designs  to  explore  the  ramifications 
of  flare  and  actuator  disk  effects.   Perhaps  it  is 
possible  to  find  a  simple  radial  equilibrium  tangential 
velocity  distribution  which  will  be  less  effected  by  the 
modification,  or  will  "compensate"  for  the  effects. 
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APPENDIX  A 

Variable  List 

The  following  is  an  alphabetical  list  of  variables 
which  appear  in  the  computer  main  program  and  subprograms 

The  units  are  the  same  as  those  listed  in  the 
Notation  section.   The  same  variable  names  are  used  to 
represent  identical  quantities  in  both  the  main  and  sub- 
programs.  In  addition,  the  subscripts  are  used  in 
consistant  manner  throughout.   An  intermediate  variable 
is  one  used  locally  in  the  program  to  simplify  arithmetic 
expressions.   Zero  is  v;ritten  0.      Unless  otherwise 
indicated,  variable  naming  follows  the  implicit  REAL/ 
INTEGER  typing  convention. 


AREA        annulus  area  between  two  adjacent 
streamlines 

A^(I)        constant  term  in  the  expression  for 
tangential  velocity 

A1(I)        coefficient  of  r  in  the  expression 
for  tangential  velocity 

A2(I)        coefficient  of  r^  in  the  expression 
for  tangential  velocity 

B(I)         coefficient  of  1/r  in  the  expression 
for  tangential  velocity 

CI  intermediate  gas  property  constant 

used  by  RNEXT 

C2  intermediate  variable 

C3  intermediate  variable 
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D 

DALPHA 
DBETA 
DD(I, J) 
DISK(I) 

DM 
DR 


DROLD 

DRZ 

DRZOLD 

F(I,1) 

F(I,2) 

F(I,3) 

F(I,4) 

FF(I,J) 

FLARE (I) 

FLOW 
FVT(I) 

FVX(I) 

G 

GAMMA 


intermediate  variable 

stator  tuning  angle 

rotor  turning  angle 

actuator  disk  correction  factor 

function  subprogram  which  calculates 
the  non-dimensional  actuator  disk 
correction  factor 

mass  flow  per  channel  betv;een 
streamlines 

difference  between  the  calculated 
and  actual  tip  radius 

value  of  DR  in  last  iteration 

difference  between  R(I,J+1)  and  RGUESS 
calculated  by  STREAM 

value  of  DRZ  in  last  iteration 
total  absolute  gas  velocity- 
total  relative  gas  velocity 
alpha;  absolute  flow  angle 
beta;  relative  flow  angle 
flare  correction  factor 

function  subprogram  which  calculates 
the  non-dimensional  flare  correction 
factor 

total  turbine  mass  flow 

function  subprogram  which  calculates 
tangential  velocity  as  a  function  of 
radius 

function  subprogram  which  calculates 
axial  velocity  as  a  simple  radial 
equilibrium  function  of  radius 

g   =  32.1740  f t-lbm/lbf /sec^ 
^o 

Y  =  ratio  of  specific  heats;  input 
constant 
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HR 
HS 

I 


ID 

I  READ 

IVJRITE 

J 

JJ 

JMl 

K 


KK 


LI 


L2 


L3 


MM 
MV 


mean  rotor  height 

mean  stator  height 

subscript  which  identifies  a  transverse 
station  along  the  turbine  axis 

1=1   stator  leading  edge 
1=2   rotor  leading  edge 
1=3   rotor  trailing  edge 

user  determined  input  data  identification 
number 

input  device  channel  number 

output  device  channel  num.ber 

streamline  index;  J  =  1,2,3,  ...,  JJ, 

starting  at  the  hub 

quantity  of  streamlines  considered 
quantity  of  flov;  channels  considered 


key  variable;   K 


1   modified  simple 
radial  equili- 
brium data 


key  variable; 


K  =  2   SRE  data 

KK  7^  1  program  diagnostics 
not  printed 

KK  =  1  diagnostics 
printed 


loop  count;  number  of  times  all  the 
strearalines  for  the  turbine  are 
calculated 

loop  count;  number  of  times  a  set  of 
streamlines  must  be  calculated  at  a 
station  to  satisfy  the  geometry 

loop  count;  the  total  number  of  times, 
during  the  calculation  of  a  set  of 
streamlines  at  a  plane,  that  the 
subprogram  RNEXT  is  called 

key  used  in  labeling  output  data 

absolute  mach  number  (REAL  variable) 
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MW 

PHIHUB 

PHITIP 

Nl 

N2 

N3 

PI 
P0(I) 

R(I,J) 

REF 

RGAS 

RGUESS 

RHO 

RNEXT(I) 

RPM 
RR 

RT(I) 
R0m 


SOLD 
SI 

TESTl 


relative  mach  number  (REAL  variable) 

hub  flare  angle 

tip  flare  angle 

test  variable  for  LI  loop  divergence 
check 

test  variable  for  L2  loop  divergence 
check 

test  variable  for  L3  loop  divergence 
check 

3.14159 

stagnation  pressure 

radial  location  of  a  streamline 

reference  radius  for  flare  correction 

specific  gas  constant;  input  data 

temporary  storage  for  a  trial  R(I,J+1) 

static  density 

function  subprogram  which  calculates 
a  trial  value  of  R(I,J+1) 

turbine  rpm;  input  data 

intermediate  variable  used  as  the 
argument  for  FLARE (I),  FVX(I),  FVT(I) 

tip  radius 

input  reference  radius  at  which  VX(I,J) 
VX0(I) 

sum  which  accumulates  the  difference 
between  streamline  positions  on  the 
current  iteration  (Ll)  and  the 
previous  (Ll-1) 

value  of  S  from  previous  iteration 

term  in  S 

test  value  to  terminate  loop  Ll 
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TEST2 
TEST3 
TEST4 

TGNT 

TR 

TJ2I(I) 

Tl 

T2 

T3 

U 

VRATIO 

VT(I,J) 

VTZ 

VX ( I , J ) 
VXR ( I ) 

VXSRE(I,J) 
VXZ 

VX0(I) 
WE 

WR 
WRATIO 

ws 
z 

Zl/  .../Z6 


test  value  to  terminate  loop  L2 

test  value  to  terminate  loop  L3 

test  value  to  ensure  the  axial 
velocity  is  positive 

tangent  of  the  flare  angle  0 

ratio  of  static  to  total  temperature 

stagnation  temperature 

input  percent  used  to  calculate  TESTl 

input  percent  used  to  calculate  TEST2 

input  percent  used  to  calculate  TESTS 

blade  speed 

tangential  velocity 

average  value  of  "tangential  velocity 
between  two  streamlines 

axial  velocity 

reference  value  of  axial  velocity; 
input  data 

simple  radial  equilibrium  axial  velocity 

average  value  of  axial  velocity  between 
two  streamlines 

input  reference  value  of  axial  velocity 

entrance  width  dimension  used  to  find 
the  flare  correction  factor 

rotor  width 

stator  v/idth 
intermediate  variable 
intermediate  variables 
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APPENDIX  B 

Program  Listing 

This  appendix  contains  the  computer  program  listing 
of  the  MAIN  program,  subroutine  STREM-l,  and  function 
subprograms  RNEXT,  FLARE,  DISK,  FVX,  and  FVT.   Refer  to 
the  table  below.   The  purpose  of  each  program  segment 
is  indicated;  flow  charts  for  MAIN,  STREAI^l,  and  RNEXT 
appear  in  Chapter  3 . 


paragraph 

segment 

page 

B.l 

MAIN 

95 

B.2 

STREAM 

103 

B.3 

functions 

106 

B.l   MAIN 

The  MAIN  program  reads  the  input  data,  governs  the 
iteration  loops  to  generate  the  streamlines,  calculates 
the  velocity  triangles,  and  prints  out  the  results.   The 
listing  begins  on  the  next  page. 
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B.2   Subroutine  STREAM 


STREAM  calculates  an  r .  .,,,  iterating  on  assumed 

1,3+1'  ^ 

values  for  V   and  V^  until  r.  .  ,  and  the  average 

X         0  1/3+1 

properties  based  on  r .  .  and  r.  .,,  are  consistent. 
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B. 3   Function  Subprograms 

The  listing  for  each  of  the  five  function  subpro- 
grams appears  on  one  of  the  following  pages. 


Function 


RNEXT 


FLARE 


DISK 


FVT 


FVX 


Page  Description 

107  given  a  value  of  r.  .  and  the 

properties  V  ,  Vg,  RNEXT 

calculates  a  trial  r.  .,, 

1/3+1 

108  given  r.  .,,,  FLARE  calculates  a 

non-dimensional  factor  which, 
multiplied  by  V   in  STREA!-!, 
gives  a  V   modified  for  FLARE 

109  given  appropriate  values  of 

V  .  . ,  DISK  calculates  a  non- 

XI, J 

dimensional  factor  used  to 
modify  axial  velocity  for 
actuator  disk  effects 

110  given  r.  .,  FVT  calculates  a  V^ 
^1/3  Q 

by  a  proposed  simple  quadratic 
relation 

111  given  r.  .,  FVX  calculates  a  V 
^1/3  ^ 

using  a  simple  radial  equili- 
brium relation 
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APPENDIX  C 

Specific  Instructions  For  Use  of  the  Program 
C.l   General 

This  appendix  is  intended  to  provide  detailed 
instructions  for  the  use  of  the  computer  program 
developed  in  this  thesis.   Preparation  of  the  input 
deck  and  interpretation  of  the  output  will  be  covered. 

The  program  was  run  on  the  Interdata  7  0  computer 
at  the  Joint  Mechanical  and  Civil  Engineering  Computer 
Facility  at  MIT,  but  should  be  easily  adaptable  to  any 
FORTRAN  IV  capable  system.   References  7  and  8  served 
as  the  FORTRAN  IV  language  guide.   All  READ  and  WRITE 
instructions  are  vi/ritten  with  dummy  variables  IREAD  and 
IWRITE,  respectively.   Therefore,  the  only  tv;o  FORTRAN 
statements  which  need  to  be  changed  in  order  to  run  on 
another  system  using  different  input/output  channel 
numbers  are  the  two  assigning  integer  values  to  IREAD 
and  IV7RITE.   The  control  cards,  listed  below,  which  may 
be  system  dependent,  are  marked  by  an  asterisk. 

C.2   Loading  Sequence 

Listed  below,  in  the  sequence  that  they  should  be 

loaded,  are  the  contituents  of  the  program: 

*  (1)   Job  card 

*(2)   Language  and  other  compilation  control 
cards 
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(3)  MAIN  program 

(4)  Subprograms  (in  any  order) 

a)  subroutine  STREAM 

b)  function  subprogram  RNEXT 

c)  function  subprogram  FLARE 

d)  function  subprogram  DISK 

e)  function  subprogram  FVT 

f)  function  subprogram  FVX 

*(5)   Execute  card 

*(6)   Common  length  designation  and  the 
execution  control  cards 

(7)   Data  cards 

*(8)   Job  termination  card 

C.3   Job  Control  Cards* 

The  program  segments  marked  by  an  asterisk  may 
vary  in  content  and  format  from  system  to  system  or  with 
time,  and  must  be  determined  from  an  appropriate  user's 
manual.   It  should  be  noted  that  the  amount  of  blank 
common  utilized  by  the  program  depends  upon  the  dimension 
of  some  of  the  variables,  which  relates  to  the  number  of 
streamlines  considered. 

C.4   Data  Cards 

The  data  deck  is  composed  of  any  number  of  separate 
data  sets  stacked  directly  together.   Each  set  is 
independent  of  every  other  set  and  represents  a  separate 
simple  radial  equilibrium  turbine  design  to  be  processed 
by  the  program. 
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Each  data  set  consists  of  ten  separate  cards  v;hich 
must  be  present  and  in  the  proper  sequence.   A  specific 
format  is  given  for  the  data  on  each  card,  although  free 
format  (data  fields  separated  by  commas)  may  be  used  if 
the  computer  system  has  that  feature. 

The  context  of  each  card  is  given  belov;  in 
paragraph  C.4.1.   A  superscript  refers  to  an  explanatory 
note  in  paragraph  C.4.2. 
C.4.1   Format  Table 

Information         Units 


Columns 

Format 

CARD 

ONE 

1-4 

14 

5-8 

14 

9-12 

14 

CARD 

Tv;o 

1-10 


FIO. 4 


11-20 

FIO.' 

CARD  THREE 

1-9 

F9.3 

10-18 

F9.3 

19-27 

F9.3 

28-36 

F9.3 

37-45 

F9.3 

46-54 

F9.3 

55-63 

F9.3 

64-72 

F9.3 

number  of  streamlines" 

diagnostics  print  key' 

3 

user  data  set  number 


gamma  (ratio  of  specific 
heats) 


none 
none 
none 

none 


specific  gas  constant     ft  lbf/lbm/°R 

turbine  rpm  rpm 

total  mass  flov;  rate       Ibm/sec 
total  temperature,  plane  1    °R 
total  pressure,  plane  1 
total  temperature,  plane  2 
total  pressure,  plane  2 
total  temperature,  plane  3 
total  pressure,  plane  3      psi 


psi 
°R 
psi 
'^R 
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CARD  FOUR 


1-10 

F10.6 

11-20 

F10.6 

21-30 

F10.6 

31-40 

F10.6 

41-50 

F10.6 

51-60 

F10.6 

CARD  FIVE 

1-10 

F10.6 

11-20 

F10.6 

21-30 

F10.6 

31-40 

F10.6 

CARD  SIX 

1-10 


11-20 


21-30 


31-40 


41-50 


51-60 


F10.6 


F10.6 


F10.6 


F10.6 


F10.6 


F10.6 


CARD  SEVEN 

1-10 

F10.5 

11-20 

F10.5 

21-30 

F10.5 

31-40 

F10.5 

hub  radius,  plane  1 
tip  radius,  plane  1 
hub  radius,  plane  3 
tip  radius,  plane  3 
stator  width 
rotor  width 

4 
test  percent  1 

test  percent  2 

test  percent  3 

test  value  4 


SRE  reference  radius, 
plane  1 

SRE  reference  velocity, 
plane  1 

SRE  reference  radius, 
plane  2 

SRE  reference  velocity, 
plane  2 

SRE  reference  radius, 
plane  3 

SRE  reference  velocity, 
plane  3 


B  for  plane  1' 


A   for  plane  1 
o     ^ 

A,  for  plane  1 
A^  for  plane  1 


ft 
ft 
ft 
ft 
ft 
ft 

none 
none 
none 
ft 

ft 

ft/sec 

ft 

ft/sec 

ft 

ft/sec 


ft  /sec 

ft/sec 

-1 
sec 

ft   sec 
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CARD  EIGHT 

(Same  as  card  seven,  but  for  plane  2) 
CARD  NINE 

(Same  as  card  seven,  but  for  plane  3) 
CARD  TEN 

(The  tenth  card  is  blank) 
C.4.2   Notes 

1.  The  limit  on  the  number  of  streamlines  is  only 
a  practical  one,  related  to  machine  computing  speed  and 
storage,  and  can  be  any  reasonable  num.ber.  The  program 
was  originally  run  with  eleven;  any  greater  number  will 
require  a  new  DIMENSION  statement  for  the  main  and  all 
subprograms.  Additional  blank  common  space  may  also 
have  to  be  allocated. 

2.  If  this  data  bit=l,  all  available  diagnostics 
will  be  printed  out.   See  paragraph  C.5.3,  below.   If  not 
equal  to  unity,  no  diagnostics  except  error  messages  will 
be  printed  out. 

3.  The  user  may  use  any  three  digit  number  to 
identify  particular  data  sets.  Some  num.ber  must  be 
supplied. 

4.  Refer  to  paragraph  3.3.13. 

5.  B,  A  ,  A, ,  A„,  are  the  coefficients  of  powers 

U       _i_       ^ 

of  r  in  the  tangential  velocity  relation.   See 
paragraph  3.3.6. 
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C.5   Output  Discription 
C.5.1   General 

The  output  for  each  data  set  may  be  considered  in 
three  parts: 

(1)  introductory 

(2)  diagnostic  (optional) 

(3)  velocity  triangle  data 

(i)   modified  simple  radial  equilibrium 
values 

(ii)  SRE  values. 

C.5. 2   Introductory  Output 

Input  data  and  other  selected  calculated  values, 
which  do  not  change  during  the  execution  of  a  data  set, 
are  printed  out  in  the  introductory  block.   The  symbols 
appearing  have  the  same  meaning  as  used  in  the  program 
and  given  in  Appendix  A.   If  the  data  is  input  also, 
the  output  units  are  the  same.   TESTl,  TEST2,  TEST3, 
REF,  and  WE  have  units  of  feet,  and  PHITIP  and  PHIHUB 
are  in  degrees. 
C.5. 3   Diagnostic  Output 

If  input  variable  KK=1,  certain  lines  of  information 
are  printed  out  during  execution  of  the  program  iterations 
Sample  output  of  this  type  is  shown  at  the  end  of  this 
appendix.   Each  line  contains  the  current  values  of 
certain  important  variables  of  the  program  at  a  certain 
step  in  the  calculations.   This  information  was  used  in 
the  development  of  the  program  and  should  not  be 
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required  routinely.   However,  it  can  increase  confidence 
in  the  results,  and  may  be  useful  if  unexpected  problems 
occur  during  operation  in  the  future.   Additional  details 
on  the  meaning  of  the  loop  counters  can  be  learned  from 
the  flow  charts  presented  in  Chapter  3. 
C.5.4   Velocity  Triangle  Output 

A  group  of  results  is  printed  for  each  of  two  cases: 

(1)  The  streamline  numbers  and  locations,  and 
velocity  results  are  printed  for  a  simple  radial  equili- 
brium preliminary  design  which  has  been  modified  by  the 
flare  and  actuator  disk  corrections  determined  by  the 
program.   Results  pertaining  to  each  of  the  planes  is 
printed,  followed  by  a  block  of  results  involving  values 
from  more  than  one  plane. 

(2)  For  comparison  purposes,  values  calculated 
at  the  same  streamlines  as  above,  but  with  no  modifica- 
tion for  flare  and  actuator  disk  are  printed  in  the 

same  format. 

The  meanings  of  the  variables  are  the  same  as 
listed  in  Appendix  A,  or  are  self-explanatory.   The  1 
and  2  appearing  in  some  colum.n  labels  refer  to  rotor 
inlet  and  outlet  triangle,  respectively,  following  the 
usual  convention.   Distances  are  in  feet,  velocities  in 
ft/sec,  and  angles  in  degrees.   The  value  VXR(I)  printed 
for  each  plane  is  the  final  value  for  the  constant  term 
in  the  FVX  relation  for  axial  velocity.   It  was  changed 
from  VX0{1)    to  satisfy  continuity. 
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C.5.5   Error  Messages 

The  input  data  to  this  program  is  assumed  to  be 
from  a  feasible  preliminary  design.   Therefore,  many 
types  of  execution  errors,  such  as  from  the  appearance 
of  negative  radii,  should  not  occur.   If  some  inconsis- 
tant  or  infeasible  data  should  be  used,  however,  measures 
have  been  taken  to  minimize  the  chances  of  hanging  the 
program  up  in  an  iteration  loop.   The  details  of  the 
testing  procedure  to  accomplish  this  are  found  in 
paragraph  3.4.   Each  iterative  loop  is  provided  with  a 
means  of  exit  if  the  error  term  increases  for  at  least 
two  successive  iterations.   Another  check  is  made  to 
determine  if  the  axial  velocity  is  less  than  1.0  ft/sec 
(an  arbitrary,  small  positive  value,  to  avoid  the 
possibility  of  dividing  by  zero) . 

If  any  of  the  conditions  listed  in  the  next  para- 
graph are  discovered,  an  error  message  will  be  printed 
out  regardless  of  whether  or  not  the  user  has  specified 
the  printing  of  the  diagnostic  data. 

The  follov/ing  is  a  list  of  error  conditions  tested 
for,  and  their  corresponding  output  messages: 

(1)   Condition:   the  axial  velocity  computed  by 
the  function  subprogram  FVX  is  less  than 
1.0  ft/sec 
Error  message:   LI  =  xxx  L2  =  xxx   L3  =  xxx 

I  =  xxx   J  =  xxx   VX(I,J)  LESS 
THAN  1,0  FT/SEC 
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ActJ-on:   execution  continues 

(2)  Condition:   the  change,  S,  in  the  1=3  plane 

streamline  positions  is  not 
decreasing  on  successive  iterations 
Error  message:   ITERATION  TERMINATED  DUE  TO 

INCREASING  VALUES  OF  S 
Action:   execution  of  current  data  set  is 
terminated 

(3)  Condition:   the  annulus  geometry  error,  DR, 

is  not  decreasing  on  successive 
iterations 
Error  message:   ITERATION  TERMINATED  DUE  TO 

INCREASING  VALUES  OF  DR 
I  =  XXX   L2  =  XXX   L3  =  xxx 
Action:   execution  of  current  data  set  is 
terminated 

(4)  Condition:   the  change,  DRZ ,  in  the  calculation 

of  R(I,J+1)  by  subroutine  STREAM 
is  not  decreasing  on  successive 
iterations 
Error  message:   ITERATION  TERMINATED  DUE  TO 

DIVERGING  POSITION  OF  STREAM- 
LINE J+1    I  =  xxx   J  =  xxx 
LI  =  xxx   L2  =  xxx   L3  =  xxx 
Action:   execution  of  current  data  set  is 
terminated 
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C.6   Additional  Options 

In  some  applications  it  may  be  desired  to  find  the 
simple  radial  equilibrium  streamlines,  or  to  consider 
the  effects  of  the  flare  or  actuator  disk  correction 
alone.   To  accomplish  this,  the  executable  statements 
in  one  or  both  of  the  subprograms  FLARE (I)  and  DISK (I) 
may  be  replaced  by  the  statement  FLARE  =  1.0  or 
DISK  =  1.0,  as  appropriate. 


121 


<■   1    .  '\  r        r.  n  z        J  J    .\j  » 

r.  ■*  (K  jf^        -t  K^  0        -  ^.i  4 

^    -»  I'  .■■         »  r-    .        r  »    /  -•, 

^•-ry  ^N.        '■ft       --  ■»-»  ^ 

^  •■•  •-•  cv  •-        ri  Ai  t\j       •->  — •  .-•  '>,. 

t  •->  iT           ■*,  J"  U     £   If     X 

r-.  ,£  (V        F  r  n  c-  r-  p. 

II    u    II   "^  M    II  -^  II    II    II  »-  n    X  II    II  i    II     .   <    X    , 

—  —  — ■r—  —  r>  —  —  —  a   —   iT—  —  r»—  *ir^*". 

"5     1    "J   ■«.  "J    T  U      "J    5     "3  ('/    -J   *  T    T  r.      ^  ^    X  r.  c> 
~~';«~-)«-lTTt-j.-3T.-;.      ... 

, ,.   »,„^_^ ^_^^f(.^ 

aaro:  .irk       naa       ti  aa  rr 

n  n  u         n  ii  n    ii    n    i> 

^►-rv  —  vr  — rr^  ^  —  .-  _r-^_-._ 

I  I                I                I  I 

^-  1."  t          fV.  (\,  #    -^   t\'   .» 

It    II    II    tik  n    «  ^   n    II    II  ■„   n    i.  II    II  ;,    u   ^  i   ti  :. 

"oo    '  c-  c  •  rz  c  c:   »  c<    '  a  Cj  •  d  '    •   •    « 

■^  re  T"i  s"^'S'5'n 

•  II  III 

o»-iinitwr"  iijr--»ii.-ii^irii*it    ii    nii 

^  ^  ^  r   -^        ^  cr  r-       \r  t  'C  ir~ 

c    r-  cv  r    r        r   rt  (v        ry  o,  fv.  n. 

r^  -?  if.        -v  c  rv    c    rv  1/ 

«    II    II    r»  u    II  X   II    II    II  r~«   «    (^  II    II  r    II  r^  ^  ii  E, 


.- 

f; 

0- 

-» 

s. 

^ 

(S. 

r^ 

^ 

r 

•c 

r 

r- 

CJ 

s 

(v 

iT. 

Ci 

CT. 

s 

X 

r^ 

r^ 

r^ 

(V 

.-• 

r 

■s 

s 

or 

M 

<* 

X 

II 

^ 

c 

-c 

II 

r' 

II 

tv 

fv 

II 

oc 

It 

» 

u 

K 

^ 

5, 

nO 

S. 

X 

s 

r^ 

ir, 

s. 

J- 

& 

in 

K 

</ 

CT. 

nC 

u: 

r- 

if- 

■x 

St 

00 

X 

. 

• 

* 

• 

. 

• 

• 

• 

• 

• 

• 

• 

hA 

^ 

p. 

•- 

-» 

*- 

-J 

a. 

p, 

r 

Jt 

r- 

^ 

o 

<£ 

<r. 

<r 

>c 

s 

•I 

a 

(^ 

0- 

a 

s 

tr 

.-. 

^^ 

w- 

M 

•-« 

Cv 

•-• 

•-• 

♦- 

»-» 

it 

»- 

a. 

•- 

cv 

jt. 

*- 

c 

ti. 

^ 

Ck 

ar 

oc 

-^ 

K 

r- 

IT 

c 

•^ 

■^ 

-' 

,"; 

11 

II 

!• 

1 

1- 

It 

0 

,, 

II 

II 

II 

„ 

a 

a 

a 

II 

a 

a 

II 

a 

a 

a 

n 

Q 

II 

a 

a 

II 

a 

II 

II 

u 

11 

„ 

> 

X 

> 

r 

> 

V 

r 

> 

X 

X 

r- 

X 

r- 

> 

X 

r 

> 

r- 

r 

T. 

r 

3   ^(vf'i-T— •fVj<^<.-itvjfo-»'-'<\j--'i\.r;«-ifV.»-'-i^ 


<     11    p    u    II    »    II    »    II    u    II    u    It    II    II    II    H    ii    II    II    II    y 
(\fv(\   fv(v<va(\rv.o.(vafvafv-aa(Vfvtv(v 

t'     _)   _l  _.    -J     J   _,     J    _i    J    _:_■_'__:_!_._.  _/    _■   _    -J 


H      U      U      II 


H     K    II    n 
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APPENDIX  D 

Integration  of  the  Simple  Radial  Equilibrium  Equation 

A  convenient  form  of  the  simple  radial  equilibrium 
equation  is 


dV     Vp,     ^ 

-V   -^  =  —  ;^(rVJ.  (D.l) 

X  dr     r  dr    e 


Given  explicit  functions,  Vg (r) ,  it  can  be  easily 
integrated  to  find  a  relation  for  axial  velocity  as  a 
function  of  radius. 

Suppose  the  expression  for  the  tangential  velocity, 
Vq  (r)  ,  can  be  written  in  the  follov^ing  form: 


Vg(r)  =  a^r^  +  a^r  +  a^  +  |.  (D.2) 


The  bottom  sign  applies  after  a  stator  row,  and  the  top 
sign  applies  after  a  rotor  row. 

Then,  preceding  to  evaluate  the  right  hand  side  of 
equation  (D. 1) , 

rVn  =  a^r^  +  a^r^  +  a  r  +  b 


A(rVQ)  =  Sa^r^  +  2a^r  +  a^ 


7^(-^e)  =    3a,r.2a,  .^;  (D.3) 
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Finally,  multiplying  (D.2)  by  (D.3),  (D.l)  becomes 
dV 


a  b 
ol    .2      'o      lr.,2 


+  (Sa^a^  +  3a_b)  +  (a  ^  +  2a, b)-  +  -2_         (d.4) 


Equation  (D.4)  may  be  integrated  directly  from  a 
reference  radius,  r.,  and  axial  velocity,  V  .,  to  an 
arbitrary  radius  r,  and  corresponding  V  .   This  gives 


+  (2a^a2  +  a^M  (r^  ~  ^i^^  "^  ^  ^^o^l  "^   ^2^^^^  "  ^i^ 


+  (a  2  +  2a,b)ln(— )  ±  a  b  (-  -  — )  . 
o      1      r .      or    r . 
1  1 


Rearranging, 


X  XI  ^2  1  ^12  1 


-    (4a   a^    +2a    2)  (r2    -   r.M    -    6  (a   a,     +   a^b)  (r    -   r. 

021  1  O     -L  ^  J- 


(2a^2    -    4a^b)ln(r/r^)     +    2a^b(|  -   p-) .  (D.5) 


The  general  expression,  (D.5),  may  be  simplified 
for  two  special  cases  in  common  use: 
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(i)  Let  a^  =  a^  =  0;  a  ,  b  7^  0.   Then  (D.5)  becomes 


V^'  -  V^^^  =  -2a^Mr'  -  r^M  ±  4a^b  ln(r/r^). 


This  represents  the  constant  reaction  type  design. 


(D.6) 
6 


(ii)  Let  a„  =  aT=0;a,b7^0.   Then  (D.5)  becomes 
2     1         o 


V  2  _  V  .2  =  -2a  2  ln(r/r.)  +  2a  b  (-  -  ~)  . 
X     XI       o        1      o   r   r. 

1 


(D.7) 


This  represents  a  design  in  which  the  stator  inlet  angle, 


a 


^,    is  approximately  the  same  at  all  radii 
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APPENDIX  E 

Actuator  Disk  Superposition 

Approximate  solutions  of  the  actuator  disk  theory, 
discussed  in  Chapter  2,    are,  in  the  case  of  the 
isolated  disk: 

(V  )    -  (V  ) 
upstream  of  the  disk  (x  <  0),  and 


(V  )    -  (V  ) 
V^  =  (V^)^„  -  -^ 2 ^—   exp(-;^)  (E.2) 


dov/nstream  of  the  disk  (x  >  0). 

These  equations  pertain  to  the  axial  variation  of 
axial  velocity  along  a  streamline. 

Suppose  it  is  desired  to  consider  the  effects  of 
these  actuator  disk  solutions  at  the  three  principal 
planes  of  a  single  stage  axial  turbine.   Two  actuator 
disks  are  involved,  corresponding  to  the  stator  row  and 
the  rotor  row.   The  three  planes  of  interest  are  the 
stator  leading  edge  (plane  1),  the  rotor  leading  edge 
(plane  2),  and  the  rotor  trailing  edge  (plane  3). 

For  simplicity,  assume  that  the  stator  leading 
edge  is  influenced  only  by  the  stator  actuator  disk 
(for  equal  blade  widths,  the  rotor  disk  is  three  times 
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the  distance  away) ,    and  that  the  rotor  trailing  edge  is 
influenced  only  by  the  rotor  actuator  disk.   The  rotor 
leading  edge  feels  the  effects  of  both  the  stator  and 
rotor  disks.   Refer  to  Figure  E-1. 

The  relations  (E.l)  and  (E.2)  apply  to  the  stage 
inlet  and  outlet  planes  respectively.   The  effect  at 
the  rotor  leading  edge,  point  S,  can  be  obtained  by 
superimposing  the  velocity  corrections  due  to  the  two 
disks.   The  effect  of  the  stator  upon  S  is  given  by 
(E.l)  v;ith  X  =  -WS/2  and  h  =  HS ;  the  effect  of  the 
rotor  is  given  by  (E.2)  with  x  =  +WR/2  and  h  =  HR. 

In  addition,  the  velocities  at  infinity,  appearing 
in  the  two  equations  must  be  interpreted.   Let  (V   ) 
denote  the  axial  velocity  a  long  way  downstream  of  the 
rotor  disk;  similarly,  let  (V   )    denote  the  axial 
velocity  a  long  way  upstream  from  the  stator,  and  so 
on.   Then  (E.l)  and  (E.2)  become 


(V   )    -  (V   ) 

V   =   V   )    +  ^ exp  -o.op  /    (E.3) 

X      xr  -oo  2  "-   2  •  HB 


which   is    the   effect   of   the   rotor,    and 


/T.         ^  ^^XS^+~  ^\s^-c«  .     TT'WS.  IT?     A\ 

V      =      V      K       -   o exp(-„^„c:)  ,  (E.4) 

X  xs    +0O  2  "■       Z*  Hb 


which  is  the  effect  of  the  stator  upon  point  S. 
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Figure  E-1 
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Now  assume 


C^vr^-co  ~  (^vc^-i-co  "^   ^^^  ~    ■the  simple  radial  equili- 
brium velocity  at  the 
rotor  leading  edge 
(plane  2)  corrected  for 
flare; 


(V   )_^  ~  ^  _  ~  "the  axial  velocity  at  plane  1; 


(V   )^   =  V  ,  =  the  axial  velocity  at  plane  3. 
xr  +«'    x+  J  ir- 


Then  finally,  the  corrected  axial  velocity,  V  ^ ,  at 
point  S  (plane  2)  can  be  written 


V         V    -  V  V    -  V 

x2  ,  ,   x+     xm     /  TT'WR,  xm     x-     ,  tt'WS. 

V—  =  1  +  — 2V ^^P^"2^Hr)  -  2V ^^P^-2^HS^ 

xm           xm  xm 


(E.5) 


V  ^  is  the  axial  velocity  corrected  for  both  flare  and 
x2  -^ 

actuator  disk  effects. 
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APPENDIX  F 

Calculation  of  the  Flare  Geometry  Factors 

The  purpose  of  this  appendix  is  to  derive  the  two 
parameters  REF  and  WE  used  by  the  computer  program  to 
compute  the  flare  angle  ^.      The  computation  of  <^    at  an 
arbitrary  point,  Q,  within  the  annulus  is  also  demon- 
strated.  The  stage  inlet  and  outlet  hub  and  shroud 
radii  are  given  information. 

Referring  to  Figure  F-1,  on  the  next  page,  it  is 
seen  that  the  point  P  is  the  intersection  of  the  extended 
line  segments  forming  the  hub  and  shroud  walls  of  the 
turbine.   The  line  x  =  0  represents  the  inlet  plane  of 
the  turbine,  x  =  x_  represents  the  outlet  plane,  and 
y  =  0  represents  the  axis  of  the  machine. 

The  equations  for  the  lines  representing  the  shroud 
and  hub  walls,  respectively,  are  v;ritten  as 


(R    -    \.) 

^■^     ^-^   X  +  R^,  (F.l) 


x^  tl 


y  =  _il3 hL-x  +  R...  (F.2) 

-^         x^  hi 


Solving  (F.l)  and  (F.2)  simultaneously  for  x  and  y,  the 
coordinates  of  the  point  P,  namely  (-WE, REF),  are  found 
to  be 
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P(-WE,REF) 


(XyR^^) 


(X3,i\3) 


>   X 


Figure  F-1 
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(R  .  -  R.  ,  )x, 
WE  =  ^li^         II       I  .  (WE  >  0)  (F.3) 

t3   tl   h3   hi 


WE 
REF  =  R,  ,  +   ^^    (R,,  -  R,^)-  (F.4) 

hi    X-    hi     h3 


The  angle  <t'  for  some  arbitrary  point,  Q,    within  the 
annulus,  can  be  found  from  the  expression 


y   -  REF 
o 


where  the  coordinates  of  the  point  Q  are  (x  ,    Y  ) • 
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